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Abstract 

A review of solar cycle prediction methods and their performance is given, including fore- 
casts for cycle 24. The review focuses on those aspects of the solar cycle prediction problem 
that have a bearing on dynamo theory. The scope of the review is further restricted to the 
issue of predicting the amplitude (and optionally the epoch) of an upcoming solar maximum 
no later than right after the start of the given cycle. 

Prediction methods form three main groups. Precursor methods rely on the value of some 
measure of solar activity or magnetism at a specified time to predict the amplitude of the 
following solar maximum. Their implicit assumption is that each numbered solar cycle is a 
consistent unit in itself, while solar activity seems to consist of a series of much less tightly 
intercorrelated individual cycles. Extrapolation methods, in contrast, are based on the premise 
that the physical process giving rise to the sunspot number record is statistically homogeneous, 
i.e., the mathematical regularities underlying its variations are the same at any point of time, 
and therefore it lends itself to analysis and forecasting by time series methods. Finally, instead 
of an analysis of observational data alone, model based predictions use physically (more or less) 
consistent dynamo models in their attempts to predict solar activity. 

In their overall performance during the course of the last few solar cycles, precursor methods 
have clearly been superior to extrapolation methods. Nevertheless, most precursor methods 
overpredicted cycle 23, while some extrapolation methods may still be worth further study. 
Model based forecasts have not yet have had a chance to prove their skills. One method that 
has yielded predictions consistently in the right range during the past few solar cycles is that 
of K. Schatten et al, whose approach is mainly based on the polar field precursor. 

The incipient cycle 24 will probably mark the end of the Modern Maximum, with the Sun 
switching to a state of less strong activity. It will therefore be an important testbed for cycle 
prediction methods, and, by inference, for our understanding of the solar dynamo. 
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1 Introduction 



Solar cycle prediction is an extremely extensive topic, covering a very wide variety of proposed 
prediction methods and prediction attempts on many different timescales, ranging from short term 
(month-year) forecasts of the runoff of the ongoing solar cycle to predictions of long term changes 
in solar activity on centennial or even millennial scales. As early as 1963, Vitinsk y published a 



whole monograph on the subject, later updat ed and extended (Vitinsk yl.ll963l . ll973fl. More recen t 



overviews of the field or aspects of it include Hathaway ( 2009f l. Kanei (|200ll ). and Pesnelll ( 2008 ). 



In order to narrow down the scope of the present review, we constrain our field of interest in two 
important respects. 

Firstly, instead of attempting to give a general review of all prediction methods suggested or 
citing all the papers with forecasts, here we will focus on those aspects of the solar cycle prediction 
problem that have a bearing on dynamo theory. We will thus discuss in more detail empirical 
methods that, independently of their success rate, have the potential of shedding some light on the 
physical mechanism underlying the solar cycle, as well as the prediction attempts based on solar 
dynamo models. 

Secondly, we will here only be concerned with the issue of predicting the amplitude (and option- 
ally the epoch) of an upcoming solar maximum no later than right after the start of the given cycle. 
This emphasis is also motivated by the present surge of interest in precisely this topic, prompted 
by the unusually long and deep recent solar minimum and by sharply conflicting forecasts for the 
maximum of the incipient solar cycle 24. 

As we will see, significant doubts arise both from the theoretical and observational side as to 
what extent such a prediction is possible at all (especially before the time of the minimum has 
become known). Nevertheless, no matter how shaky their theoretical and empirical backgrounds 
may be, forecasts must be attempted. Making verifiable or falsifiable predictions is obviously the 
core of the scientific method in general; but there is also a more imperative urge in the case of 
solar cycle prediction. Being the prime determinant of space weather, solar activity clearly has 
enormous technical, scientific, and financial impact on activities ranging from space exploration 
to civil aviation and everyday communication. Political and economic decision makers expect the 
solar community to provide them with forecasts on which feasibility and profitability calculations 
can be based. Acknowledging this need, the Space Weather Prediction Center of the US National 
Weather Service does present annually or semiannually updated "official" predictions of the up- 
coming sunspot maximum, emit ted by a Solar Cycle Prediction Panel of experts, starting shortly 
before the (expected) minimum (jSWPCt ). The unusual lack of consensus during the early meetings 



of this panel during the recent minimum, as well as the con current more frequently updated but 
wildly varying predictions of a NASA MSFC team (MSFC) have put on display the deficiencies 



of currently applied prediction techniques; on the other hand, they also imply that cycle 24 may 
provide us with crucial new insight into the physical mechanisms underlying cyclic solar activity. 

While a number of indicators of solar activity exist, by far the most commonly employed is still 
the smoothed relative sunspot number R; the "Holy Grail" of sunspot cycle prediction attempts 
is to get R right for the next maximum. We therefore start by briefly introducing the sunspot 
number and inspecting its known record. Then, in Sections [21 El and U we discuss the most widely 
employed methods of cycle predictions. Section [5] presents a summary evaluation of the past 
performance of different forecasting methods and collects some forecasts for cycle 24 derived by 
various approaches. Finally, Section [6] concludes the paper. 

1.1 The sunspot number 

Despite its somewhat arbitrary construction, the series of relative sunspot numbers constitutes 
the longest homogeneous global indicator of solar activity determined by direct solar observations 
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and carefully controlled methods. For this reason their use is still predominant in studies of solar 
activity variation. As defined originally by IWoll ( 1850l ) . the relative sunspot number is 



R w = k(10g + f) (1) 

where g is the number of sunspot groups (including solitary spots), / is the total number of all spots 
visible on the solar disc, while k is a correction factor depending on a variety of circumstances, such 
as instrument parameters, observatory location and details of the counting method. Wolf, who 
decided to count each spot only once and not to count the smallest spots, the visibility of which 
depended on seeing, used k — 1. The counting system employed was changed by Wolf's successors 
to count even the smallest spots, attributing a higher weight (i.e. / > 1) to spots with a penumbra, 
depending on their size and umbral structure. As the new counting naturally resulted in higher 
values, the correction factor was set to k = 0.6 for subsequent determinations of Rw to ensure 
continuity with Wolf's work, even though there was no change in either the instrument or the 
observing si t e. Th is w as followed by several f urther changes in the deta ils of the counting meth od 
( Waldmeier . 1961 ; see Kopeckv et al , 1980l . Hovt and Schatten . 19981 and Hathaway . feoiObl for 



further discussions on the determination of Rw)- 

In addition to introducing the relative sunspot number. IWoll (1861) also used earlier observa- 



tional records available to him to reconstruct its monthly mean values since 1749. In this way, 
he reconstructed 11-year sunspot cycles back to that date, introducing their still universally used 
numbering. (In a later work he also determined annual mean values for each calendar year going 
back to 1700.) 

In 1981, the observatory responsible for the official determination of the sunspot number 
changed from Zurich to the Royal Observatory of Belgium in Brussels. The website of the SIDC 
(originally Sunspot Index Data Center, recently renamed Solar Influences Data Analysis Center), 
idc . oma . be , is now the most authoritative source of archive sunspot number data. But it 
has to be kept in mind that the sunspot number is also regularly determined by other institutions: 
these variants are informally known as the American sunspot number (collected by AAVSO and 
available from the National Geophysical Data Center, http://www.ngdc.noaa.gov/ngdc.html) 
and the Kislovodsk Sunspot Number (available from the web page of the Pulkovo Observatory, 



http://www.gao.spb.ru). Cycle amplitudes determined by these other centers may differ by up 



to 6-7% from the SIDC values, NOAA numbers being consistently lower, while Kislovodsk numbers 
show no such systematic trend. 

These significant disagreements between determinations of Rw by various observatories and 
observers are even more pronounced in the case of historical data, especially prior to the mid-19th 
century. In particular, the controversial suggestion that a whole solar cycle may have been missed 
in the official sunspot number series at the end of the 18th century is taken by some as glaring 
evidence for the unreliability of early observations. Note, however, that independently of whether 
the claim for a missing cycle is well founded or not, there is clear evidence that this controversy 
is mostly due to the very atypical behaviour of the Sun itself in the given period of time, rather 
than to the low quality and coverage of contemporary observations. These issues will be discussed 
further in section [3.2.21 

Given that Rw is subject to large fluctuations on a time scale of days to months, it has become 
customary to use annual mean values for the study of longer term activity changes. To get rid of 
the arbitrariness of calendar years, the standard practice is to use 13-month boxcar averages of the 
monthly averaged sunspot numbers, wherein the first and last months are given half the weight of 
other months: 




R — — + 2 Rm,i + Rm,6 J (2) 

R m ,i being the mean monthly value of Rw for ith calendar month counted from the present month. 
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It is this running mean R that we will simply call "the sunspot number" throughout this review 
and what forms the basis of most discussions of solar cycle variations and their predictions. 

(n) (n) 

In what follows, -Rmax and R m [ n will refer to the maximum and minimum value of R in cycle n 
(the minimum being the one that starts the cycle). Similarly, imax and tj™ n will denote the epochs 
when R takes these extrema. 



1.1.1 Alternating series and nonlinear transforms 

Instead of the "raw" sunspot number series R(t) many researchers prefer to base their studies on 
some transformed index R' . The motivation behind this is twofold. 

(a) The strongly peaked and asymmetrical sunspot cycle profiles strongly deviate from a si- 
nusoidal profile; also the statistical distribution of sunspot numbers is strongly at odds with a 
Gaussian distribution. This can constitute a problem as many common methods of data analy- 
sis rely on the assumption of an approximately normal distribution of errors or nearly sinusoidal 
profiles of spectral components. So transformations of R (and, optionally, t) that reduce these 
deviations can obviously be helpful during the analysis. In this vein, e.g., Max Waldmeier often 
based his studies of the solar cycle on the use of logarithmic sunspot numbers R' = log i?; many 
other researchers use R' = R a with 0.5 < a < 1, the most common value being a = 0.5. 

(b) As the sunspot number is a rather arbitrary construct, there may be an underlying more 
physical parameter related to it in some nonlinear fashion, such as the toroidal magnetic field 
strength B, or the magnetic energy, proportional to B 2 . It should be emphasized that, contrary to 
some claims, our current understanding of the solar dynamo does not make it possible to guess what 
the underlying parameter is, with any reasonable degree of certainty. In particular, the often used 
assumption that it is the magnetic energy, lacks any sound foundation. If anything, on the basis 
of our current best understanding of flux emergence we might expect that the amount of toroidal 
flux emerging from the tachocline should be J \B — Bq \ dA where Bq is some minimal threshold 
field strength for Par ker instabili t y and the surface integral goes across a latitudinal cross section 
of the tachocline (cf. iRuzmaikin , 1997f). As, however, the lifetime of any given sunspot group i s 



finite and proportional to its size ( Petrovav and van Driel-Gesztelvi . 1997 ; Henwood et ai , 2010h . 



instantaneous values of R or the total sunspot area should also depend on details of the probability 
distribution function of B in the tachocline. This just serves to illustrate the difficulty of identifying 
a single physical governing parameter behind R. 

One transformation that may still be well motivated from the physical point of view is to 
attribute an alternating sign to even and odd Schwabe cycles: this results in the the alternating 
sunspot number series R±. The idea is based on Hale's well known polarity rules, implying that 
the period of the solar cycle is actually 22 years rather than 11 years, the polarity of magnetic 
fields changi ng sign from one 11-year Schwabe cycle to the next. In this representation, first 
suggested by iBracewell (1953), usually odd cycles are attributed a negative sign. This leads to 



slight jumps at the minima of the Schwabe cycle, as a consequence of the fact that for a 1-2 year 
period around the minimum, spots belonging to both cycles are present, so the value of R never 
reaches zero; in certain applications, further twists are introduced into the transformation to avoid 
this phenomenon. 

After first introducing the alternating series, in a later work IBracewell (Il988l) demonstrated 
that introducing an underlying "physical" variable Rb such that 

R ± = 100(i? B /83) 3/2 (3) 

(i.e. a = 2/3 in the power law mentioned in in item (a) above) significantly simplifies the cycle 
profile. Indeed, upon introducing a "rectified" phase variably 4> in each cycle to compensate for 

1 The more precise condition defining <f> is that <j> = ±n/2 at each maximum and <f> is quadratically related to the 
time since the last minimum. 
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the asymmetry of the cycle profile, Rb is a nearly sinusoidal function of <j). The empirically found 
3/2 law is interpreted as the relation between the time-integrated area of a typical sunspot group 
vs. its peak area (or peak Rw value), i.e. the steeper than linear growth of R with the underlying 
physical parameter Rb would be due to to the larger sunspot groups being observed longer, and 
therefore giving a disproportionately larger contribution to the annual mean sunspot numbers. If 
this interpretation is correct, as suggested by Bracewell's analysis, then Rb should be considered 
proportional to the total toroidal magnetic flux emerging into the photosphere in a given interval. 
(But the possibility must be kept in mind that that the same toroidal flux bundle may emerge 
repeatedly or at different heliographic longitudes, giving rise to several active regions.) 



1.2 Other indicators of solar activity 



Reconstructions of R prior to the early 19th century are increasing ly uncertain. In order to 
tackle problems related to sporadic and often unreliable observations, iHovt and Schatten (1998) 
introduced the Group Sunspot Number (GSN) as an alternative indicator of solar activity. In 
contrast to Rw, the GSN only relies on counts of sunspot groups as a more robust indicator, 
disregarding the number of spots in each group. Furthermore, while Rw is determined for any 
given day from a single observer's measurements (a hierarchy of secondary observers is defined for 
the case if data from the primary observer were unavailable), the GSN uses a weighted average 
of all observations available for a given day. The GSN series has been reproduced for the whole 
period 1611-1998 (Figure [I} and it is generally agreed that for the period 1611-1818 it is a more 
reliable reconstruction of solar activity than the relative sunspot number. Yet there have been 
relatively few attempts to date to use this data series for solar cycle prediction. One factor in this 
could be the lack of regular updates of the GSN series, i.e. the unavailability of precise GSN values 
for the past decade. 
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Figure 1: 13- mo nth sliding averages of the monthly average relative sunspot numbers R (green) 
and group sunspot numbers Rq (black) for the period 1611-1998. 

Instead of the sunspot number, the total area of all spots observed on the solar disk might 
seem to be a less arbitrary measure of solar activity. However, these data have been available since 
1874 only, covering a much shorter period of time than the sunspot number data. In addition, 
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the determination of sunspot areas, especially farther from disk center, is not as trivial as it may 
seem, resulting in significant random and systematic errors in the area determinations. Area 
measurements performed in two different observatories often sh ow discrepancies reach ing ~ 30% 
for smaller spots (cf. the figure and discussion in Appendix A of Petrovav et al. , 19991 ). 



A number of other direct indicators of solar activity have become available from the 20th 
century. These include e.g., various plage indices or the 10.7 cm solar radio flux - the latter is 
considered a particularly good and simple to measure indicator of global activity (cf. Figure [2J ■ 
As, however, these data sets only cover a few solar cycles, their impact on solar cycle prediction 
has been minimal. 




1950 1960 1970 1980 1990 2000 

Year 



Figure 2: Monthly values of the 10.7 cm radio flux in solar flux units for the period 1947-2009. 
The solar flux unit is defined as 10~ 22 W/m 2 Hz. The green curve shows R m + 60 where R m is the 
monthly mean relative sunspot number. (The vertical shift is for better comparison.) Data are 
from the NRC Canada (Ottawa/Penticton). 



Of more importance are proxy indicators such as geomagnetic indices (the most widely used 
of which is the aa index), the occurrence frequency of aurorae or the abundances of cosmogcnic 
radionuclides such as 14 C and 10 Be. For solar cycle prediction uses such data sets need to have 
a sufficiently high temporal resolution to reflect individual 11-year cycles. For the geomagnetic 
indices such data have been avail able since 1868. while an annual 10 Be series covering 600 years has 
been published very recently by Berggren et al. ( 20091 ). Attempts have been made to reconstruct 



the epochs and even amplitudes of solar maxima d uring the past two millennia from oriental naked 
eye s unspot records and from auroral observations ([Stephenson and Wolfendald . Il988l ; iNagovitsvnl . 



1997), but these reconstructions are currently subject to too many uncertainties to serve as a basis 
for predictions. Isotopic data with lower temporal resolution are now available for up to 50,000 
years in the past; while such data do not show individual Schwabe cycles, they are still useful for 
the study of long term variations in cycle amplitude. Inferring solar activity parameters from such 
proxy data is generally not straightforward. 
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Figure 3: The variation of the monthly smoothed relative sunspot number R during the period 
1749-2009, with the conventional numbering of solar cycles. 
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1.3 The solar cycle and its variation 

The series of R values determined as described in Section 11.11 is plotted in figure [3J It is evident 
that the sunspot cycle is rather irregular. The mean length of a cycle (defined as lasting from 
minimum to minimum) is 11.02 years (median 10.7 years), with a standard deviation of 1.2 years. 
The mean amplitude is 113 (median 115), with a standard deviation of 40. It is this wide variation 
that makes the prediction of the next cycle maximum such an interesting and vexing issue. 

It should be noted that the period covered by the relative sunspot number record includes 
an extended interval of atypically strong activity, the so called Modern Maximum (see below), 
cycles 17-23. Excluding these cycles from the averaging, the mean and median values of the cycle 
amplitude are very close to 100, with a standard deviation of 35. The mean and median cycle 
length then become 11.1 and 11.2 years, respectively, with a standard deviation of 1.3 years. 



1.3.1 Secular activity variations 



Inspecting Figure [3] one can discern an obvious long term variation. For the study of such long 
term variations, the series of cycle parameters is often smoothed on time scales significantly longer 
than a solar cycle: this procedure is kn own as secular sm oothing. One popular method is the 
so-called Gleissberg filter or 12221 filter ( Gleissberel . 1967 ). For instance, the Gleissberg filtered 
amplitude of cycle n is given by 



(«) 
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Figure 4: Amplitudes of the sunspot cycles (dotted) and their Gleissberg filtered values (blue 
solid), plotted against cycle number. 

The Gleissberg filtered sunspot number series is plotted in Fig. 4. One long-term trend is an 
overall secular increase of solar activity, the last six or seven cycles being unusually strong. (Four 
of them are markedly stronger than average and none is weaker than average.) This period of 
elevated sunspot activity level from the mid-20th century is known as the "Modern Maximum". 
On the other hand, cycles 5, 6 and 7 are unusually weak, forming the so-called "Dalton Minimum". 
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Finally, the rather long series of moderately weak cycles 12-16 is occasionally referred to as the 
"Gleissberg Minimum" - but note that most of these cycles are less than la below the long-term 
average. 

While the Dalton and Gleissberg minima are but local minima in the ever changing Gleissberg 
filtered SSN series, the conspicuous lack of sunspots in the period 1640-1705, known as the Maunder 
Minimum (Figure [T]) quite obviously represents a qualitatively different state of solar activity. 
Such extended periods of high and low activity are usually referred to as grand maxima and 
grand minima. Clearly, in comparison with the Maunder Minimum, the Dalton Minimum could 
only be called a "semi-grand minimum" , while for the Gleissberg Minimum even that adjective is 
undeserved. 

A number of possibilities have been proposed to explain the phe nomenon of grand minima and 



maxima, including chaotic behaviou r of the non l inear solar dynamo (IWeiss et all . Il984l ). stochastic 



fluctuations in dynamo parameters (IMoss et all 120081 ; lUsoskin et a i. 2009b ) or a bimodal dynamo 



with stochastically induced alternation between two stationary states (jPetrovavl . 120071 ). 

The analysis of long-term proxy data, extending over several millennia further showed that 
there exist systematic long-term statistical trends and periods such as the so called secular and 
supersecular cycles (see Section EO]) . 



1.3.2 Does the Sun have a long term memory? 

Following customary usage, by "memory" we will refer to some physical (or, in the case of a model, 
mathematical) mechanism by which the state of a system at a given time will depend on its previous 
states. In any system there may be several different such mechanisms at work simultaneously — if 
this is so, again following common usage we will speak of different "types" of memory. A very 
mundane example are the RAM and the hard disk in a computer: devices that store information 
over very different time scales and the effect of which manifests itself differently in the functioning 
of the system. 

There is no question that the solar dynamo (i.e. the mechanism that gives rise to the sunspot 
number series) does possess a memory that extends at least over the course of a single solar cycle. 
Obviously, during the rise phase solar activity "remembers" that it should keep growing, while in 
the decay phase it keeps decaying, even though exactly the same range of R values are observed in 
both phases. Furthermore, profiles of individual sunspot cyc les may, in a first approximation, be 
considered a one-parameter ensemble ( Hathaway et aZ.I . Il994h . This obvious effect will be referred 



to here as intracycle memory. 

As we will see, correlations between activity parameters in different cycles are generally much 
weaker than those within one cycle, which strongly suggests that the intracycle memory mechanism 
is different from longer term memory effects, if such are present at all. Referring back to our 
analogy, the intracycle memory may work like computer RAM, periodically erased at every reboot 
(i.e. at the start of a new cycle). 

The interesting question is whether, in addition to the intracycle memory effect, any other type 
of memory is present in the solar dynamo or not. To what extent is the amplitude of a sunspot 
cycle determined by previous cycles? Are subsequent cycles essentially independent, randomly 
drawn from some stochastic distribution of cycle amplitudes around the long term average? Or, 
in the alternative case, for how many previous cycles do we need to consider solar activity for 
successful forecasts? 

The existence of long lasting grand minima and maxima suggests that the sunspot number 
record must have a long-term memory extending over several consecutive cycles. Indeed, elemen- 
tary combinatorical calculations show that the occurrence of phenomena like the Dalton minimum 
(3 of the 4 lowest maxima occurring in a row) or the Modern maximum (4 of the 5 highest max- 
ima occurring within a series of 5 cycles) in a random series of 24 recorded solar maxima has a 
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rather low probability (5% and 3%, respectively). This conclusion is corroborated by the analysis 
of long-term proxy data, extending over several millennia, which showed that the occurrence of 
grand minima and gr and maxima is more common than what would follow from Gaussian statistics 
(jUsoskin et ZI . I2007I ). 



It could be objected that for sustained grand minima or maxima a memory extending only from 
one cycle to the next would suffice. In contrast to long-term (multidecadal or longer) memory, this 
would constitute another kind of short-term (Ss 10 years) memory: a cycle-to-cycle or intercycle 
memory effect. In our computer analogy, think of system files or memory cache written on the hard 
disk, often with the explicit goal of recalling the system status (e.g. desktop arrangement) after 
the next reboot. While these files survive the reboot, they are subject to erasing and rewriting in 
every session, so they have a much more temporary character than the generic data files stored on 
the disk. 

The intercycle memory explanation of persistent secular activity minima and maxima, however, 
would imply a good correlation between the amplitudes of subsequent cycles, which is not the case 
(cf. Sect. 12.11 below). With the known poor cycle-to-cycle correlation, strong deviations from the 
long-term mean would be expected to be damped on time scales short compared to e.g. the length 
of the Maunder minimum. This suggests that the persistent states of low or high activity are due 
to truly long term memory effects extending over several cycles. 

Further evidence for a long-term memory in solar activity comes from the persistence analysis 
of activity indicators. The parameter determined in such studies is the Hurst exponent < H < 1. 
Essentially, H is the steepness of the growth of the total range 1Z of measured values plotted 
against the number n of data in a time series, on a logarithmic plot: 1Z oc n H . For a Markovian 
random process with no memory H ~ 0.5. Processes with H > 0.5 are persistent (they tend to 
stay in a stronger-than-average or weaker-than- average state longer), while those with H < 0.5 are 
anti-persistent (their fluctuations will change sign more often). 

Hurst expo nents for solar activity indices have been derived using rescaled range analysis by 
many authors dMandelbrot and Wallisl. Il969l: iRuzmaikin et aH Il994l: Komm. 1995; Oliver and 



Ballester, 1 19961 : iKilcik et all 120091 ). All studies coherently yield a value H = 0.85-0.88 for time 
scales exceeding a year or so, and somewhat lower values (H ~ 0.75) on shorter time scales. Some 
doubts regarding the significance of this result for a finite scries have been raised bv Oliver and 



Ballester (|1998f ): however. lOian and Rasheedl (|2004 ) have shown using Monte-Carlo experiments 



that for time series of a length comparable to the sunspot record, H values exceeding 0.7 are 
statistically significant. 

A complementary method, essentially equi valent to rescaled range analysis is detrended fluc- 



tuation analysis. Its application to solar data (jOgurtsovi . 120041 ) has yielded results in accordance 
with the H values quoted above. 

The overwhelming evidence for the persistent character of solar activity and for the intermittent 
appearance of secular cyclicities, however, is not much help when it comes to cycle-to-cycle predic- 
tion. It is certainly reassuring to know that forecasting is not a completely idle enterprise (which 
would be the case for a purely Markovian process) , and the long-term persistence and trends may 
make our predictions statistically somewhat different from just the long-term average. There are, 
however, large decadal scale fluctuations superposed on the long term trends, so the associated 
errors will still be so large as to make the forecast of little use for individual cycles. 

1.3.3 Waldmeier effect and amplitude frequency correlation 



Greater activity on the Sun goes with shorter periods, and less with longer 
periods. I believe this law to be one of the most important relations among the 
Solar actions yet discovered. 

(jWoli Il86lh 
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Figure 5: Monthly smoothed sunspot number R at cycle maximum plotted against the rise time 
to maximum (left) and against cycle length (right). Cycles are labeled with their numbers. In the 
plots the red dashed lines are linear regressions to all the data, while the blue solid lines are fits to 
all data except outliers. Cycle 19 is considered an outlier on both plots, cycle 4 on the right hand 
plot only. The corresponding correlation coefficients are shown. 



It is apparent from Figure [3] that the profile of sunspot cycles is asymmetrical, the rise being 
steeper than the decay. Solar activity maxima occur 3 to 4 years after the minimum, while it 
takes another 7-8 years to reach the next minimum. It can also be noticed that the degree of this 
asymmetry correlates with the amplitude of the cycle: to be more specific, the length of the rise 
phase anticorrelates with the maximal value of R (Figure [5]), while the length of the decay phase 
shows weak or no such correlation. 

Historically, the relation was first formulated by IWaldmeierl f|l935h as an inverse correlation 
between the rise time and the cycle amplitude; however, as shown by iTritakisl (1982), the total 
rise time is a weak (inverse logarithmic) function of the rise rate, so this representation makes 
the correlation appear less robust. (Indeed, when formulate d with the rise time it is not even 
present in some activity indicators, suc h as sunspot areas - cf. Dikpati et al. (2008b).) As pointed 
out bv ICameron and Schussler ( 20081 ). the weak link between rise time and slope is due to the 
fact that in steeper rising cycles the minimum will occur earlier, thus partially compensating for 
the shortening due to a higher rise rate. T he effec t is indeed more clearly seen when th e rate 
of the rise is used instead of the rise time (|Lantosl . [2000l : 1(5 ameron and Schiisslerl . 120081 ). The 
observed correlation between rise rate and maximum cycle amplitude is approximately linear, 
good (correlation coefficient r ~ 0.85), and quite robust, being present in various activity indices. 

Nevertheless, when coupled with the nearly nonexistent correlation between the decay time and 
the cycle amplitude, even the weaker link between the rise time and the maximum amplitude is 
sufficient to forge a weak inverse correlation between the t otal cycle le ngth and the cycle amplitude 
(Figure [5]). This inverse relationship was first noticed by IWolJ(ll86ll) . 

A stronger i nverse correlation was found between the cycle amplitude and the length of the 
previous cycle by lHathaway et ~al. (1993). This correlation is also r eadily explained as a cons e quenc e 
of the Waldmeier effect, as demonstrated in a simple model bv ICameron and Schiisslerl (|2007l ). 
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Note that in a more detailed study ISolanki et all (|2002l ) find that the correlation coefficient of this 
relationship has steadily decreased during the course of the historical sunspot number record, while 
the correlation between cycle amplitude and the length of the third preceding cycle has steadily 
increased. The physical significance (if any) of this latter r es ult is unclear. 

I n what follows, the relationships found bv Wolfi (Il86lh . lHathawav et all ( 19941) . and Solanki 
et al. (|2002l ). discussed above, will be referred to as "R^^-try^ r correlations" with n = 0. —1 or 
—3, respectively. 

Modern time series analysis methods offer several ways to define an instantaneous frequency / in 
a quasiperiodic series . One simple approac h was discussed in the context of Bracewell's transform, 
Equation ([3]), above. iMininni et al. ( 2000f ) discuss several more sophisticated methods to do this, 
concluding that Gabor's analytic signal approach yields the best p erformance. This technique was 
first applied to the sunspot record bv IPalus and Novotnal (1999), who found a significant long 
term correlation between the smoothed instantaneous frequency and amplitude of the signal. On 
time scales shorter than the cycle length, however, the frequency-amplitude correlation has not 
been convincingly proven, and the fact that the correlation coefficient is close to the one reported 
in the right hand panel of Figure [5] indicates that all that the fashionable gadgetry of nonlinear 
dynamics could achieve was to recover the effect already known to Wolf. It is clear from this that 
the "frequency-amplitude correlation" is but a secondary consequence of the Waldmeier effect. 

On the left hand panel of Figure [5j within the band of correlation the points seem to be 
sitting neatly on two parallel strings. Any number of faint hearted researchers wo uld dism i ss this 
as a coincidence or as another manifestation of the "Martian canal effect". But iKuklinl (jl98fih 
boldly speculated that the phenomenon may be real. Fair enough, cycles 22 and 23 dutifully took 
their place on the lower string eve n after the publica tion of Kuklin's work. This speculation was 
supported with further evidence bv iNagovitsvn (1997) who offered a physical explanation in terms 
of the amplitude-frequency diagram of a forced nonlinear oscillator (cf. Section 14.51) . 

Indeed, an anticorrelation between cycle length and amplitude is characteristic of a class of 
stochastically forced nonlin ear oscillators and it may also be reproduced by introducing a stochasti c 
forcing in dynamo models (Stixj, Il972l : lOssendriiver et 12, Il996l : ICharbonneau and Dikpatl l2000h . 
In some such model s the characteristic asymmetric profile of the cycle is also well reproduced 
( Mininni et Mi l2000l l2002h . The predicted amplitude-frequency relation has the form 



logR&l = C X + C 2 f. 



(5) 



Nonlinear dynamo models including some form of a-quenching also have the potential to re- 
produce the effects described by Wolf and Waldmeier wit hout recourse to stochastic driving . In 
a dynamo with a Kleeorin-Ruzmaikin type feedback on a, iKitiashvili and Kosovichev (2009) arc 
able to qualitatively reproduce the Waldmeier effect. Assuming that the sunspot number is related 
to the toroidal field strength according to the Bracewell transform, Equation ([3]) , they find a strong 
link between rise time and amplitude, while the correlations with fall time and cycle length are 
much weaker, just as the observations suggest. They also find that the form of the growth time- 
amplitude relationship differs in the regular (multiperiodic) and chaotic regimes. In the regular 
regime the plotted relationship suggests 



R {n) =Ci -Co( & L) - t (n) ) 

"max u l I I ""max t min J ' 



while in the chaotic case 



'/(' 



I/ I An) _M) 
1 1 "max mm 



Note that based on the actual sunspot number series Waldmeier originally proposed 



io gJ RL n 2x = Ci-c 2 (41x-il) 



(6) 
(7) 

(8) 
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while according to Dmitrieva et al. ( 20001 ) the relation takes the form 



]QttB&l<x 



(n) 



(9) 



At first glance, these logarithmic empirical relationships seem to be more compatible with 
the relation (J5J) predicted by the stochastic models. These, on the other hand, do not actually 
reproduce the Waldmeier effect, just a general asymmetric profile and an amplitude-frequency 
correlation. At the same time, inspection of the the left hand panel in Figure [5] shows that the 
data is actually not incompatible with a linear or inverse rise time-amplitude relation, especially if 
the anomalous cycle 19 is ignored as an outlier. (Indeed, a logarithmic representation is found not 
to improve the correlation coefficient - its only advantage is that cycle 19 ceases to be an outlier.) 
All this indicates that nonlinear dynamo models may have the potential to provide a satisfactory 
quantitative explanation of the Waldmeier effect, but more extensive comparisons will need to be 
done, using various models and various representations of the relation. 
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2 Precursor Methods 



Jeder Fleckenzyklus muB als ein abgeschlossenes Games, als ein Phiinomen fur 
sich, aufgefaBt werden, und es reiht sich einfach Zyklus an Zyklus. 



dGleissberd . Il952h 



In the most general sense, precursor methods rely on the value of some measure of solar ac- 
tivity or magnetism at a specified time to predict the amplitude of the following solar maximum. 
The precursor may be any proxy of solar activity or other indicator of solar and interplanetary 
magnetism. Specifically, the precursor may also be the value of the sunspot number at a given 
time. 

In principle, precursors might also herald the activity level at other phases of the sunspot cycle, 
in particular the minimum. Yet the fact that practically all the good precursors found need to be 
evaluated at around the time of the minimum and refer to the next maximum is not simply due to 
the obvious greater interest in predicting maxima than predicting minima. Correlations between 
minimum parameters and previous values of sol ar indices have been looked for, but the results 
were overwhelmingly negative (e.g. Tlatov . 2009). This indicates that the sunspot number series 



is not homogeneous and Rudolf Wolf's instinctive choice to start new cycles with the minimum 
rather than the maximum in his numbering system is not arbitrary - for which even more obvious 
evidence is provided by the butterfly diagram. Each numbered solar cycle is a consistent unit in 
itself, while solar activity seems to consist of a series of much less tightly intercorrelated individual 
cycles, as suggested by Wolfgang Gleissberg in the motto of this section. 

In Section 11.3.21 we have seen that there is significant evidence for a long-term memory under- 
lying solar activity. In addition to the evidence reviewed there, systematic long-term statistical 
trends and periods of solar activity, such as the secular and supersecular cycles (to be discussed in 
Section I3.2|) also attest to a secular mechanism underlying solar activity variations and ensuring 
some degree of long-term coherence in activity indicators. However, as we noted, this long-term 
memory is of limited importance for cycle prediction due to the large, apparently haphazard decadal 
variations superimposed on it. What the precursor methods promise is just to find a system in 
those haphazard decadal variations — which clearly implies a different type of memory. As we 
already mentioned in Section ll.3.2[ there is obvious evidence for an intracycle memory operating 
within a single cycle, so that forecasting of activity in an ongoing cycle is currently a much more 
successful enterprise than cycle-to-cycle forecasting. As we will see, this intracycle memory is one 
candidate mechanism upon which precursor techniques may be founded, via the Waldmeier effect. 

The controversial issue is whether, in addition to the intracycle memory, there is also an intercy- 
cle memory at work, i.e. whether behind the apparent stochasticity of the cycle-to-cycle variations 
there is some predictable pattern, whether some imprint of these variations is somehow inherited 
from one cycle to the next, or individual cycles are essentially independent. The latter is known 
as the "outburst hypothesis" : consecutive cycles would then represent a series of "outbursts" of 
activity with stochastica lly fluctuating amplitudes (jHalml . Il90ll : IWaldmeieJ . Il935l : IVitinskyl . 1973; 



see also Ide Meverl . Il98ll who calls this "impulse model" ) . Note that cycle-to-cycle predictions in 
the strict temporal sense may be possible even in the outburst solar cycles are known to 

overlap. Active regions belonging to the old and new cycles may coexist for up to three years or so 
around sunspot minima; and high latitude ephemeral active regions orient ed according to the next 
cycle appear as early as 2-3 years after the maximum ( Tlatov et al , 201dl - the so-called extended 
solar cycle). 

In any case, it is undeniable that for cycle-to-cycle predictions, which are our main concern 
here, the precursor approach seems to have been the relatively most successful, so its inherent 
basic assumption must contain an element of truth - whether its predictive skill is due to a "real" 
cycle-to-cycle memory (intercycle memory) or just to the overlap effect (intracycle memory). 
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The two precursor types that have received most attention are polar field precursors and geo- 
magnetic precursors. A link between these two categories is forged by a third group, characterizing 
the interplanetary magnetic field strength or "open flux" . But before considering these approaches, 
we start by discussing the most obvious precursor type: the level of solar activity at some epoch 
before the next maximum. 




Figure 6: Monthly smoothed sunspot number R at cycle maximum plotted against the values of 
R at the previous minimum (left) and 2.5 years before the minimum (right). Cycles are labeled 
with their numbers. The blue solid line is a linear regression to the data; corresponding correlation 
coefficients are shown. In the left hand panel, cycle 19 was considered an outlier. 



2.1 Cycle parameters as precursors and the Waldmeier effect 

The simplest weather forecast method is saying that "tomorrow the weather will be just like 
today" (works in about 2/3 of the cases). Similarly, a simple approach of sunspot cycle prediction 
is correlating the amplitudes of consecutive cycles. There is indeed a marginal correlation, but 
the correlation coefficient is quite low (0.35). The existence of the correlation is related to secular 
variations in solar activity, while its weakness is due to the significant cycle-to-cycle variations. 

A significantly better correlation exists between the minimu m activity lev el and the the am- 
plitude of the next maximum (Figure [5]) . The relation is li near (|Brownl. 1976 ) , with a c orrelation 
coeffi cient of 0.72 (if the anomalous cycle 19 is ignored - Braisa et al , 20091 : see also Pishkalo . 
20081 ). The best fit is 

R max = 67.5 + 6.91i? min (10) 

Using the observed value 1.7 for the SSN in the recent minimum, the next maximum is predicted 
by this "m inimax " method to rea ch values around 80, with a la error of about ±25. 

Cameron and Schussler (|2007l ) point out that the activity level three years before the minimum 
is an even better predictor of the next maximum. Indeed, playing with the value of time shift we 
find that the best correlation coefficient corresponds to a time shift of 2.5 years, as shown in. the 
right hand panel of Figure [S] (but this may depend on the particular time period considered, so we 



15 



will refer to this method in Table 1 as "minimax3" for brevity) . The linear regression is 

i?max = 41.9 + lMR{t min - 2.5). (11) 

For cycle 24 the value of the predictor is 16.3, so this indicates an amplitude of 69, suggesting that 
the upcoming cycle may be comparable in strength to those during the Gleissberg minimum at the 
turn of the 19th and 20th centuries. 

As the epoch of the minimum of R cannot be known with certainty until about a year after 
the minimum, the practical use of these methods is rather limited: a prediction will only become 
available 2-3 years before the maximum, and even then with the rather low reliability reflected in 
the correlatio n coe fficients qu oted above. In additio n, as convincing ly demonst rated by Cameron 
and Schiissler (|2007l ) in a Monte Carlo simulation, these methods do not constitute real cycle-to- 
cycle prediction in the physical sense: instead, they are due to a combination of the overlap of solar 
cycles with the Waldmeier effect. As stronger cycles are characterized by a steeper rise phase, the 
minimum before such cycles will take place earlier, when the activity from the previous cycle has 
not yet reached very low levels. 

The same overlap readily explains the i? m ax _ icycie,n correlations discu ssed in Sect ion 11.3.31 
These relationships may also be used for solar cycle prediction purposes (e.g. Kanel. 2008h but they 
lack robustness. For cycle 24 the -R ma x-icycie,-i correlation, as formu lated bv iHathawayl (|2010tJ ) 
predicts i? ma x = 80 while the methods used bv ISolanki et al. ( 20021 ) yield values ranging from 
86 to about 110, depending on the relative weights of i Gy cie,-i an d i C ycie,-3- The forecast is not 
only sensitive to the value of n used but also to the data set (relative or group sunspot numbers) 
( Vaquero and Trigol . Ejoih . 
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Figure 7: Magnetic field strength in the Sun's polar regions as a function of time. Blue solid: North; 
red dashed: ( — 1) -South; thin black solid: average; heavy black solid: smoothed average). Strong 
annual modulations in the hemispheric data are due to the tilt of the solar equator to the Ecliptic. 



Data and figure: Wilcox Solar Observatory - see htt p : //wso . Stanford . edu/gif s/Polar . gif | for 
updated version 
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2.2 Polar precursors 

Direct measurements of the magnetic field in the polar areas of the Sun hav e been available 
from Wilcox Observatory since 1976 (pvalgaard et aL , llQTSl: iHoeksemal . Il995h . Even before a 



significant amount of data had been available for statistic al analysis, sole l y on the basis of the 
Babcock-Leighton scenario of the origin of the solar cycle, ISchatten et al. ( 1978 ) suggested that 



the polar field measurements may be used to predict the amplitude of the next solar cycle. Data 
collected in the four subsequent solar cycles have indeed confirmed this suggestion. As it was 
originally motivated by theoretical considerations, this polar field precursor method might also be 
a considered a model-based prediction technique. As, however, no particular detailed mathematical 
model is underlying the method, numerical predictions must still be based on empirical correlations 
- hence our categorization of this technique as a precursor method. 

The shortness of the available direct measurement series represents a difficulty when it comes to 
finding empirical correlations to solar activ ity. This problem can to some extent be circumvented 
by the use of proxy data. For instance, lObridko and Shelting (2008) use H Q synoptic maps 



to reconstruct the polar field strength at the source surface back to 1915. Spherical harmonic 
expansions of global photospheric magnetic measurements can also be used to deduce the field 
strength near the poles. The use of such proxy techniques permits a forecast with a sufficiently 
restricted error bar to be made, despite the shortness of the direct polar field data set. 

The polar fields reach their maximal amplitude near minima of the sunspot cycle. In its most 
commonly used form, the polar field precursor method employs the value of the polar magnetic 
field strength (typically, the absolute value of the mean field strength poleward of 55° latitudes, 
averaged for the two hemispheres) at the time of sunspot minimum. It is indeed remarkable that 
despite the very limited available experience, forecasts using the polar field method ha ve proven 
to b e consistently in the right range for cycles 21, 22 and 23 (ISchatten and Sofia . 19871; Schatten 



et aL ll996l ). 

By virtue of the definition (|2|), the time of the minimum of R cannot be known earlier than 
6 months after the minimum - indeed, to make sure that the perceived minimum was more than 
just a local dip in R, at least a year or so needs to elapse. This would suggest that the predictive 
value of polar field measurements is limited, the prediction becoming available 2-3 years before 
the upcoming maximum only. 

To remedy this situation, ISchatten and Pesnelll (1993) introduced a new activity index, the 



"Solar Dynamo Amplitude" (SoDA) index, combining the polar field strength with a traditional 
activity indicator (the 10.7 cm radio flux F10.7). Around minimum, SoDA is basically proportional 
to the polar precursor and its value yields the prediction for i<T0.7 at the next maximum; however, 
it was constructed so that its 11-year modulation is minimized, so theoretically it should be rather 
stable, making predictions possible well before the minimum. That's the theory, anyway - in reality, 
SoDA based forecasts made more than 2-3 years before the minimum usually proved unreliable. 
It is then questionable to what extent SoDA improves the prediction skill of the polar precursor, 
to which it is more or less equivalent in those late phases of the solar cycle when forecasts start to 
become reliable. 

Fortunately, however, the maxima of the polar field curves are often rather flat (cf. Figure [7]), 
so approximate forecasts are feasible several years befo re the actual min i mum . Using the current, 
rather flat and low maximum in polar field strength, ISvalgaard et al. ( 2005 ) have been able to 



predict a relatively weak cycle 24 (peak R value 75 ± 8) as early as 4 years before the sunspot 
minimum took place in December 2008! Such an early prediction is not always possible: early polar 
field predictions of cycles 22 and 23 had to be corrected later and only forecasts made shortly before 
the actual minimum did finally converge. Nevertheless, even the moderate success rate of such 
early predictions seems to indicate that the suggested physical link between the precursor and the 
cycle amplitude is real. 
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In addition to their above mentioned use in reconstructing the polar field strength, various 
proxies or alternative indices of the global solar magnetic field during the activity minimum may 
also be used directly as activity cycle precursors. H Q synoptic charts are now available from various 
observatories from as early as 1870. As H a filaments lie on the magnetic neutral lines, these maps 
can be used to reconstruct t he overall topology, if not the detailed map, of the large-scale solar 
magnetic field. iTlatovl (|2009h has shown that several indices of the polar magnetic field during the 
activity minimum, determined from these charts, correlate well with the amplitude of the incipient 
cycle. 

High resolution Hinode observations have now demonstrated that the polar magnetic field has 
a strongly int ermittent structure, b eing concentrated in intense unipolar tubes that coincide with 
polar faculae ( Tsuneta et al . 20081 ). The number of polar faculae should then also be a plausible 
proxy of the polar magnetic field strength and a good pre cursor of t he inc ipient solar cycle around 
the minimum . This conclusion was indeed confirmed by iLi et al. (2002), and, more recently, by 
Tlatovi (|2009l 'l. 

These methods offer a prediction over a time span of 3-4 years, comparable to the rise time of 
the next cvcle. A significantly earlier prediction possibilitv was. however, suggested bv Makarov 



et al. (|l989l ) and iMakarov and Makaroval (|l996l ) based on the number of polar faculae observed 
at Kislovodsk, which was found to predict the next sunspot cycle with a time lag of 5-6 years; 
even short term annual variations or "surges" of sunspot activity were claimed to be discernible in 
the polar facular record. This surprising result may be partly due to the fact that Makarov et al. 
considered all faculae poleward of 50 degrees latitude. Bona fide polar faculae, seen on Hinode 
images to be knots of the unipolar field around the poles, are limited to higher latitudes, so the 
wider sample may consist of a mix of such "real" polar faculae and small bipolar ephemeral active 
regions. These latte r are known to obey an extended butterfly diagram, as recently confirmed by 
Tlatov et al] (|201dh : the first bipoles of the new cycle appear at higher latitudes about 4 years 
after the activity maximum. It is not i mpossible that these ea rly ephemeral active regions may 
be used for predictio n purposes (cf. also Badalvan et all . 2001); but whethe r or not the result of 
Mak arov et al. ljl989h may be attributed to this is doubtful, as Ei et all (|2002h find that even using 
all polar faculae poleward of 50° from the Mitaka data base, the best autocorrelation still results 
with a time shift of about 4 years only. 

Finally, trying to correlate various solar activity parameters, Tlatov ( 20091 ) finds this surprising 
relation: 



tf££ X) =Ci + C 2 R^l (t£l - 4"2x) d = 83 ± 11 C 2 = 0.09 ± 0.02 



(12) 



where irev is the epoch of the polarity reversal in cycle n (typically, about a year after t 



(n) 



The 



origin of this curious relationship is unclear. In any case, the good correlation coefficient (r = 0.86, 
based on 12 cycles) and the time lag of ~ 10 years make this relationship quite remarkable. For 



cycle 24 this formula predicts R 



(24) 



94 ±14 



2.3 Geomagnetic and interplanetary precursors 

Relations between the cycle related variations of geomagnetic indices and solar activity were noted 
long ago. It is, however, important to realize that the overall correlation between geomagnetic 
indices and solar activity, even after 13-month smoothing, is generally far from perfect. This is 
due to the fact that the Sun can generate geomagnetic disturbances in two ways: 

(a) By material ejections (such as CMEs or flare particles) hitting the terrestrial magnetosphcrc. 
This effect is obviously well correlated with solar activity, with no time delay, so this con- 
tribution to geomagnetic disturbances peaks near, or a few years after, sunspot maximum. 
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(Note that the occurrence of the largest flares and CMEs is known to peak some years after 



the sunspot maximum - see Figure 16 in lHathawavl . 12010b 



(b) By a variation of the strength of the general interplanetary magnetic field and of solar wind 
speed. Geomagnetic disturbances may be triggered by the alternation of the Earth's crossing 
of interplanetary sector boundaries (slow solar wind regime) and its crossing of high speed 
solar wind streams while well within a sector. The amplitude of such disturbances will clearly 
be higher for stronger magnetic fields. The overall strength of the interplanetary magnetic 
field, in turn, depends mainly on the total flux present in coronal holes, as calculated from 
potential field source surface models of the coronal magnetic field. At times of low solar 
activity the dominant contribution to this flux comes from the two extended polar coronal 
holes, hence, in a simplistic formulation this interplanetary contribution may be considered 
linked to the polar magnetic fields of the Sun, which in turn is a plausible precursor candidate 
as we have seen in the previous subsection. As the polar field reverses shortly after sunspot 
maximum, this second contribution often introduces a characteristic secondary minimum in 
the cycle variation of geomagnetic indices, somewhere around the maximum of the curve. 

The component (a) of the geomagnetic variations actually follows sunspot activity with a 
variable time delay. Thus a geomagnetic precursor based on features of the cycle dominated by 
this component has relatively little practical u tility. This would seem to be the case, e.g., with the 



forecast method first proposed by l0h.ll (|1966l ). who noticed that the minimum amplitudes of the 
smoothed geom agnetic aa index are correlated to the amplitude of the next sunspot cycle (see also 
Du et allmM ). 



An indication that the total geomagnetic activity, resulting from both mecha nisms does contai n 
useful information on the expected amplitude of the next solar cycle was given bv lThompson (1993), 



who found that the total number of disturbed days in the ge omagnetic field in cycle n is related 
to the sum of the amplitudes of cycles n and n + 1. (See also lDabas et al. . 20081 ) . 



A method for separating component (b) was proposed bv lFevnman (|l982h who correlated the 



annual aa index with the annual mean sunspot number and found a linear relationship between R 
and the minimal value of aa for years with such R values. She interpreted this linear relationship 
as representing the component (a) discussed above, while the amount by which aa in a given year 
actually exceeds the value predicted by the linear relation would be the contribution of type (b) 
(the "interplanetary component"). The interplanetary component usually peaks well ahead of 
the sunspot minimum and the amplitude of the peak seemed to be a good predictor of the next 
sunspot maximum. However, it is to be noted that the assumption that the "surplus" contribution 
to aa originates from the interplanetary component only is likely to be erroneous, especially for 
stronger cycles. It is known that the number of large solar eruptions shows no unique relation to 
R: in p articular, for R > 100 their frequency may vary by a factor of 3 (see Figure 15 in lHathawavl . 
2010bl ). so in some years they may well yield a contribution to aa that greatly exceeds the minimum 



contribution. A case in point was the "Halloween events" of 20 03, that ve r y likely resulted in a large 
false contribution to the derived "interplanetary" aa index (|Hathawavl . 2010d). As a result, the 



geomagnetic precursor method based on the separation of the interplanetary component predicts 
an unusually strong cycle 24 (R m ~ 150), in contrast to most other methods, including Ohl's 
method and the polar field precursor, which suggest a weaker than average cycle (R m ~ 80-90). 

In addition to the problem of neatly separating the interplanetary contribution to geomagnetic 
disturbances, it is also wrong to assume that this interplanetary c ontribution is dominated by the 
effect of polar magnetic fields at all times during the cycle. Indeed, IWang and Sheelevl (|2009l ) point 



out that the interplanetary magnetic field amplitude at the Ecliptic is related to the equatorial 
dipole moment of the Sun that does not survive into the next cycle, so despite its more limited 
practical use, Ohl's original method, based on the minima of the aa index is physically better 
founded, as the polar dipole dominates around the minimum. The total amount of open interplan- 
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etary flux, more closely linked to polar fields, could still be determined from geomagnetic activity 
if the interplanetary contribution to it is further split into 

(bl) A contribution due to the varying solar wind speed (or to the interplanetary magnetic field 
strength anticorrelated with it), which in turn reflects the strength of the equatorial dipole. 

(b2) Another contribution due to the overall interplanetary field strength or open magnetic flux, 
which ultimately reflects the axial dipole. 

Clearly, if the solar wind speed contribution (bl) could also subtracted, a physically better founded 
prediction method should result. While in situ spacecraft measurements for the solar wind speed 
and the in terp lanetar y magnetic field strength do not have the necessary time coverage. Svalgaard 
and Oliver (|2005l . l2007h and lRouillard et all (l2007f ) devised a method to reconstruct the variations of 



both v ariables from geomagnetic measurements alone. Building on their results. I Wang and Sheelev 



(2009) arrive at a prediction of R m = 97 ± 25 for the maximum amplitude of solar cycle 24. To 



what extent the effect of the Halloween 2003 events has been removed from this analysis is unclear. 
In any case, the prediction agrees fairly well with that of Bhatt et al. (2009) who, assuming a 



preliminary minimum time of August 2008 and applying a modified form of Ohl's method, predict 
a cycle maximum in late 2012, with an amplitude of 93 ± 20. 

The actual run of cycle 24 will be certainly most revealing from the point of view of these 
complex interrelationships. 

The open magnetic flux can also be derived from the extrapolation of solar magnetograms using 
a potential field source surface model. The magnetograms applied for this purpose may be actual 
observations or the output from surface flux transport models, using the sunspot distribution 
(butterfly diagram) and the meridional flow as input. Such models indicate that the observed 
latitude independence of the interplanetary field strength ("split monopole" structure) is only 
reproduced if the source surface is far enough (> 10 Rq) and the potential field model is modified 
to take into a ccount the he liospher ic current sheet (current sheet sourc e surface model, Schiissler 
and Baumann, 20061 : Jiang et al . 2010al ). The extrapolations are generally found to agree well 



with in situ measurements where these are available. 
2.4 Flows in the photosphere 

In the currently widely popular flux transport dynamo models the strong polar fields prevalent 
around sunspot minimum are formed by the advection of following polarity flux from active regions 
by the poleward meridional flow. Changes in this flow may thus influence the would-be polar fields 
and thereby may serve as precursors of the upcoming cycle. 

Such changes, on the other hand, are also associated with the normal course of the solar activity 
cycle, the overall flow at mid-latitudes being slower before and during maxima and faster during 
the decay phase. Therefore, it is just the cycle-to-cycle variation in this normal pattern that may 
be associated with the activity variations between cycles. In this respect it is of interest to note 
that the poleward flow in the late phases of cycle 2 3 seems to have had an excess speed relative to 



the previous cycle (jHathawav and Rightmird . 120101 ) . If this were a latitude- independent amplitude 



modulation of the flow, then most flux transport dynamo models would predict a stronger than 
average polar field at t he minimum, contr ary to observations. On the other hand, in the surface 
flux transport model of Wang et al. (|2009l ) an increased poleward flow actually results in weaker 



polar fields, as it lets less leading polarity flux to d iffuse across the equator and cancel there. As 
the recent analysis bv lMunoz-Jaramillo et al'. ( 2010l ) has shown, the discrepancy resulted from the 



neglect of leading polarity flux in the Babcock-Leighton source term in flux transport dynamo 
models, and it can be remedied by substituting a pair of opposite polarity flux rings as source term 
instead of the a-term. With this correction, 2D flux transport and surface flux transport models 
agree in predicting a weaker polar field for faster meridional flow. 
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It is known from helioseismology that meridional flow speed fluctuations follow a characteristic 
latitudinal pattern associated with torsional oscillations and the butterfly diagram, consisting of 
a pair of axisymmetric bands of latitudinal flows converging towards the activity belts, migrating 
towards the equator, and accompanied by similar high-latitude poleward branches. This suggests 
interpreting the unusual meridional flow speeds observed during cycle 23 as an increased amplitude 
of this migrating modulation, rather than a change in the large-scale flow speed (Cameron and 
Schussler. l2010l ). In this case, the flows converging on the activity belts tend to inhibit the transport 
of following polariti es to the poles, agai n resulting in a lowe r than usual polar field, as observed 
( Jiang et all [2010b ; note, however, that Svanda et al. , 2007 find no change in the flux transport 
in areas with increased flows). It is interesting to note that the torsional oscillation pattern, and 
thus presumably the associated meridional flow modulation pattern, was shown to be fairly well 
reproduced by a microquench i ng m echanism due to magnetic flux emerging in the active belts 
( Petrovav and Forgacs-Daikal . [20021 ). Observational support for this notion has been provided 



by th e seismic detection of locally increased flow modulation near active regions (| Svanda et al 



20071 ). This suggests that stronger cycles may be associated with a stronger modula tion pattern 



introdu cing a nonlinearity into the flux transport dynamo model, as suggested by Ijiang et al 
( 2010bl ). 

In addition to a variation in t he amplitude of m igrating flow modulations, their migration 
speed may also influence the cycle. iHowe et al. ( 20091 ) point out that in the current minimum the 
equatorward drift of the torsional oscillation shear belt corresponding to the active latitude of the 
cycle has been slower than in the previous minimum. They suggest that this slowing may explain 
the belated start of cycle 24. 
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3 Extrapolation Methods 



In contrast to precursor methods, extrapolation methods only use the time series of sunspot num- 
bers (or whichever solar activity indicator is considered) but they generally rely on more than one 
previous point to identify trends that can be used to extrapolate the data into the future. They 
are therefore also known as time series analysis, or, for historic reasons, regression methods. 

A cornerstone of time series analysis is the assumption that the time series is homogeneous, 
i.e., the mathematical regularities underlying its variations are the same at any point of time. 
This implies that a forecast for, say, three years ahead has equal chance of success in the rising or 
decaying phase of the sunspot cycle, across the maximum, or, in particular, across the minimum. 
In this case, distinguishing intracycle and intercycle memory effects, as we did in Sections 11.3. 21 and 
[2J would be meaningless. This concept of solar activity variations as a continuous process stands 
in contrast to that underlying precursor methods, where solar cycles are thought of as individual 
units lasting essentially from minimum to minimum, correlations within a cycle being considerably 
stronger than from one cycle to the next. While, as we have seen, there is significant empirical 
evidence for the latter view, the possibility of time homogeneity cannot be discarded out of hand. 
Firstly, if we consider the time series of global parameters (e.g., amplitudes) of cycles, homogeneity 
may indeed be assumed fairly safely. This approach has rarely been used for the directly observed 
solar cycles as their number is probably too low for meaningful inferences - but the long data sets 
from cosmogenic radionuclides are excellent candidates for time series analysis. 

In addition, there may be good reasons to consider the option of homogeneity of solar activity 
data even on the scale of the solar cycle. Indeed, in dynamo models the solar magnetic field simply 
oscillates between (weak) poloidal and (strong) toroidal configuration: there is nothing inherently 
special about either of the two, i.e., there is no a priori reason to attribute a special significance 
to solar minimum. While at first glance the butterfly diagram suggests that starting a new cycle 
at the minimum is the only meaningful way to do it, there may be equally good arguments for 
starting a new cycle at the time of polar reversal. There is, therefore, plenty of motivation to try 
and apply standard methods of time series analysis to sunspot data. 

Indeed, as the sunspot number series is a uniquely homogeneous and long data set, collected 
over centuries and generated in a fairly carefully controlled manner, it has become a favorite testbed 
of time ser ies analysis methods and is routinely used in t extbooks and monographs for illustration 



purposes. ([Box and Jenkins! . 120081 : IWeil . 120051 : iTond . Il993h . This section will summarize the various 



approaches, proceeding, by and large, from the simplest towards the most complex. 
3.1 Linear regression 

Linear (auto)regression means representing the value of a time series at time t by a linear combi- 
nation of values at times t — At, t — 2At, . . ., t — pAt. Admitting some random error e„, the value 
of R in point n is 

p 

Rn = Ro + J]] c„_i_R„_i + e„ 

i=l 

where p is the order of the autoregression and the Cj's are weight parameters. A further twist on 
the model admits a propagation of errors from the previous q points: 

p i 

Rn -^0 ~t~ ^ ^ C-n—iRn—i ~t~ ^ ^ ^n — i^n— i 

i=l i=l 

This known as the ARMA (AutoRegressive Moving Average) model. 
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Linear regression techniques have been widely used for solar activity prediction during the 
course of an o ngoing cycle. Their application for cycle-to-cycle predict i on ha s been less common 



course ot an ongoi ng cycle, ineir application tor cycle-to-cy 
and successful (TLpmb and Anderse nt ll980l :lB ox and Jenkins! 



20081 : Fv^ . l2005h 



Brajsa et al. ( 20091 ) applied an ARMA model to the series of annual values of R. A successful 
fit was found for p = 6, q = 6. Using this fit, the next solar maximum was predicted to take place 
around 2012.0 with an amplitude 90 ± 27, and the following minimum occurring in 2017. 



ng 

Instead of applying an autoregression model directly to SSN data, iHiremathl (|2008l ) applied it 
to a forced and damped harmonic oscillator model claimed to well represent the SSN series. This 
resulted in a predicted amplitude of 110 ± 10 for solar cycle 24, with the cycle starting in mid-2008 
and lasting 9.34 years. 



3.2 Spectral methods 



...the use of any mathematical algorithm to derive hidden periodicities from 
the data always entails the question as to whether the resulting cycles are 
not introduced either by the particular numerical method used or by the time 
interval analyzed. 

ride Meverl . Il98lh 



Spectral analysis of the sunspot number record is used for prediction under the assumption 
that the main reason of variability in the solar cycle is a long-term modulation due to one or more 
periods. 

The usual approach to the problem is the purely formal one of representing the sunspot record 
with the superposition of eigenfunctions forming an orthogonal basis. From a technical point of 
view, spectral methods are a complicated form of linear regression. The analysis can be performed 
by any of the widely used means of harmonic analysis: 

(1) Least squares (LS) frequency analysis (sometimes called "Lomb-Scargle periodogram" ) 
consists in finding by trial and error the best fitting sine curve to the data using the least squares 
method, subtracting it ( "prewhitening" ) , then repeating the procedure until the residuals become 
indistinguishab le from white noise. The first serious attempt at sunspot cycle prediction, due to 
Kimural (|l913T ). belonged to this group. The analysis resulted in a large number of peaks with 
dubious physical significance. The prediction given for the upcoming cycle 15 failed, the forecasted 
amplitude being ~ 60 while the cycle actually peaked at 105. However, it is interesting to note that 
Kimura correctly predicted the long term strengthening of solar activity duri ng the first half of the 
20th century! LS frequency analysis on sunspot data was also performed bv lLomb and Andersen 
(1980), with similar results for the spectrum. 

(2) Fourier analysis is probably the most commonly used method of spectral decom position 
in scie nc e. It has been ap p lied to sunspot data from the beginning of the 20th century ([Turnerl . 
1913al lbl: iMichelsonl . [l9ll . IVitinskvl (|l973h judg es Fourier-b ased forecasts even less reliable than 
LS periodogram methods. Indeed, for instance ICold (1973) predicted cycle 21 to have a peak 
amplitude of 60, while the real value proved to be nearly twice that. 

(3) The maximum entropy method (MEM) relies on the Wiener-Khinchin theorem that the 
power spectrum is the Fourier transform of the autocorrelation function. Calculating the autocor- 
relation of a time series for M <^ N points and extrapolating it further in time in a particular way 
to ensure maximal entropy can yield a spectrum that extends to arbitrarily low frequencies despite 
the shortness of the data segment considered, and also has the property of being able to reproduce 
sharp spectral features (if such are present in the data in the first place). A good d escription of 



the m ethod is given by I Abies (1974), accompanied with some propaganda for it - see Press et al 



(1992) for a more balanced account of its pros and cons. The use of MEM for sunspot number 
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prediction was pioneered by ICurriel (1973). Using maximum entropy me thod combin ed with mul- 
tiple regression analysis (MRA) to estimate the amplitudes and phases, iKand (|2007h arrived at a 
prediction of 80 to 101 for th e max i mum amplitude of cycle 24. It should be noted that the same 
method yielded a prediction (|Kanel . Il999f ) for cycle 23 that was far off the mark. 

(4) Singular spectrum analysis ( SSA ) is a relatively novel method for the orthogonal decompo- 
sition of a time series. While in the methods discussed above the base was fixed (the trigonometric 
functions), SSA allows for the identification of a set of othogonal eigenfunctions that are most 
suitable for the problem. This is done by a principal componen t analysis of t he co variance matrix 
Tik = {RiRi+k}- SSA was first applied to the sunspot recor d by Rangaraianl (ll998T) who only used 
this method for pre-filtering before the application of MEM. lLoskutov et al\ (|200ll ) who also give a 
good description of the method, already made a prediction for cycle 24: a peak amplitude of 117 . 



More recently, the forecast has been corrected slightly downwards to 106 ( Kuzanvan et al . 2008) 



The dismal performance of spectral predictions with the methods (l)-(3) indicates that the 
sunpot number series cannot be well represented by the superposition of a limited number of fixed 
periodic components. Instead, 

• the periods may be time dependent 

• the system may be quasiperiodic, with a significant finite width of the periodic peaks (esp. 
the 11-year peak) 

• there may be non-periodic (i.e., chaotic or stochastic) components in the behaviour of the 
system, manifested as a continuous background in the spectrum. 

In practice, all three effects suggested above may play some part. The first mentioned effect, 
time dependence, may in fact be studied within the framework of spectral analysis. MEM and 
SSA are intrinsically capable of detecting or representing time dependence in the spectrum, while 
LS and Fourier analysis can study time dependence by sliding an appropriate data window across 
the period covered by observations. If the window is Gaussian with a width proportional to the 
frequency we arri ve at the popula r wave le t analysis. This method was applied to the s unspo t 



number series bv lOchadlick et al. (Il993l) . Vigouroux and Delachie ( 19941 ). Frick et al. ( 1997 ). 



Fligge et al] (jl999j) and 



Li et all 



2005) who could co nfirm the existence and s light variation of 



the 11-year cycle and the Gleissberg-cycle. Recently, iKollath and Olahl (|2009h called attention 
to a variety of other generalized time dependent spectral analysis methods, of which the pseudo- 
Wigner transform yields especially clear details (cf. Figure [9]) . The time varying character of the 
basic periods makes it difficult to use these results for prediction purposes but they are able to 
shed some light on the variation as well as the presistent or intermittent nature of the periods 
determining solar activity. 

In summary, it is fair to say that forecasts based on harmonic analysis are notoriously unreliable. 
The secular variation of the basic periods, obeying as yet unknown rules, would render harmonic 
analysis practically useless for the prediction of solar cycles even if solar activity could indeed be 
described by a superposition of periodic functions. Although they may be potentially useful for 
very long term prediction (on centennial scales), when it comes to cycle-to-cycle forecasts the best 
we can hope from spectral studies is apparently an indirect contribution, by constraining dynamo 
models with the inambiguously detected periodicities. 

In what remains from this subsection, we briefly review what these apparently physically real 
periods are and what impact they may have on solar cycle prediction. 



3.2.1 The 11-year cycle and its harmonics 

As an example of the period spectrum obtained by these methods, in Figure [5] we present the FFT 
based power spectrum estimate of the smoothed sunspot number record. Three main features are 
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Figure 8: Power spectrum of the smoothed monthly sunspot number series for the period 1749- 
2008. Solid vertical bars mark the 11-year period, its first harmonic and subharmonic; dashed 
vertical bars are drawn at a fiducial period of 14.5 years, its harmonic and subharmonic. 



immediately noticed: 

• The dominant 11-year peak, with its sidelobes and its 5.5-year harmonic 

• The 22-year subharmonic, representing the even-odd rule 

• The significant power present at periods longer than 50 years, associated with the Gleissberg 
cycle 

The dominant peak in the power spectrum is at ~ 11 years. Significant power is also present 
at the first harmonic of this period, at 5.5 years. This is hardly surprising as the sunspot number 
cycles, as presented in Figure EH have a markedly asymmetrical profile. It is a characteristic of 
Fourier decomposition that in any periodic series of cycles where the profiles of individual cycles 
are non-sinusoidal, all harmonics of the base period will appear in the spectrum. 

Indeed, were it not for the 13-month smoothing, hig her harmonics could a l so be expected to 
appear in the power spectrum. It has been proposed ( Krivova and Solankil . 2002 ) that these 
harmonics are detected in the sunspot record and that t hey may be related to the periodicities of 



~ 1.3 years intermittent l y observed in solar wind sp eed ( Richardso n et all Il994t iPaularena et 
19951 ; ISzabo et all ll995l; iMursula and Ziegel l200ol ; [f~ 



Lockwood 



_____ . 111 ) and in the internal rotation 

velocity of the Sun (iHowd. I2009T Sect. 10.1). An analoguous intermittent 2.5 year variation in 
the solar neutrino flux (jShirail . .20041 ) may also belong to this group of phenomena. It may be 
worth noting tha t, from the o t her end of the period spect rum, the 154-day Rieger period in solar 
flare occurrence ( Rieger et all Il984 ; Bai and Oliver . 1990l ) has also been tentatively linked to the 
1.3-year periodicity. Unusually strong excitation of such high harmonics of the Schwabe cycle may 
possi bly be ex plained by excitation due to unstab le Ross bv waves in the ta chocline (Zaqarashvili 
et _L l2010l .. 

The 11-year peak in the power spectrum has substantial width, related to the rather wide 
variation in cycle lengths in the range 9-13 years. Yet Figure [5] seems to suggest the presence of a 
well detached second peak in the spectrum at a period of ~ 14 years. The presence of a distinct 
peak at the first harmonic and even at the subharmonic of this period seems to support its reality. 
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Indeed, peaks at a round 14 and 7 years were already found by other researchers (e.g., iKimura 



19131 : ICurrid . 119731 ) who suggested that these may be real secondary periods of sunspot activity. 

The situation is, however, more prosaic. Constraining the time interval considered to data more 
recent than 1850, from which time the sunspot number series is considered to be more reliable, 
the 14.5-year secondary peak and its harmonics completely disappear. On the other hand, the 
power spectrum for the years 1783-1835 indicates that the appearance of the 14.5-year secondary 
peak in the complete series is almost entirely due to the strong predominance of this period (and 
its harmonic) in that interval. This interval covers the unusually long cycle 4 and the Dalton 
minimum, consisting of three consecutive unusually weak cycles, when the "normal" 11-year mode 
of operation was com pletely suppress ed. 



>pcration was completely suppressed. 

As pointed out bv lPetrovay ( 2010h . this probably does not imply that the Sun was operating in 



a different mode during the Dalton minimum, the cycle length being 14.5 years instead of the usual 
11 years. Instead, the effect may be explained by the well known inverse correlation between cycle 
length and amplitude, which in turn is the consequence of the strong inverse correlation between 
rise rate and cycle amplitude (Waldmeier effect), combined with a much weaker or nonexistent 
correlation between decay rate and amplitude (see Section fl .3 .3[) . The cycles around the Dalton 
minimum, then, seem to lie at the low amplitude (or long period) end of a continuum representing 
the well known cycle length-amplitude relation, ultimately explained by the Waldmeier effect. 

A major consequence of this is that the detailed distribution of peaks v aries significantly de- 
pending on the interval of time considered. Indeed, Kollath and Olahl ( 20091) recently applied time 



dependent harmonic analysis to the sunspot number series and found that the dominant periods 
have shown systematic secular changes during the past 300 years (Figure [9|). For instance, the 
basic period seems to have shortened from 11 years to 10 years between 1850 and 1950, with some 
moderate increase in the last 50 years. (This is consistent with the known anticorrelation between 
cycle length and amplitude, cf. Section [1.3.31 ) 

3.2.2 The even odd (a.k.a. Gnevyshev Ohl) rule 

A cursory look at Figure [3] shows that solar cycles often follow an alternating pattern o f higher 



and lo wer maxima. In this apparent pattern, already noticed by the early observers (e.g, I Turner 



1913d ). odd cycles have been typically stronger than even cycles in the last two centuries. 

This even-odd rule can be given two interpretations: a "weak" one of a general tendency of 
alternation between even and odd cycles in amplitude, or a "strong" one of a specific numerical 
relation between the amplitudes of consecutive cycles. 

Let us first consider the rule in its weak interpretation. At first sight the rule admits many 
exceptions, but the amplitude of solar cycles depends on the particular measuring method used. 
Exceptions from the even-odd alternation rule become less common if a long term tr end (calculated 
by ap plying a 12221 or 121 filter, see Section fl . 3 . 1 [) is subtracted from the data ( Charbonneaul . 
20011 ) , or if integrated cycle amplitu des (sums of annual mean sunspot numbers during the cycle) 



are used (jGnevvshev and Ohlll94. 



In fact, as evident from, e.g., the work of lMursula et al\ (|200l[ ) where cycle amplitudes are 



based on group sunspot numbers and the amplitude of a cycle is defined as the sum of the annual 
GSN value over the course of the cycle, the odd-even alternation may be considered as strictly 
valid with only four exceptions: 

• In the pairs 7-8 and 17-18, odd cycles are followed by stronger even cycles at the end of 
Dalton minimum and at the beginning of the Modern Maximum. These excep tions could be 
made to disappear by the subtraction of the long term trend as suggested by Charbonneaul 
(|200lh . 

• The pair 22-23 represents another apparent break of the weak even-odd rule which is not 
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Figure 9: Pseudo-Wigner power distribution in the sunspot number record, with time on the 
abscissa and frequency on the ordinate. The three horizontal bands of high power correspond, 
from bottom to top, to the Gleissberg cycle, the 11-year cycle and its first harmonic. The sunspot 
number curve is shown on top for guidance. (Figure courtesy Z. Kollath) 
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easily explained away, even though the re l ative difference is smaller if the Kislovodsk sunspot 



number series is used ( Nagovitsvn et al . 20091 ). The possibility is obviously there that the 



subtraction of the long term trend may resolve the problem but we have no way to tell in 
the near future. 

• Prior to cycle 5, the phase of the alternation was opposite, even cycles being stronger than 
odd cycles. As cycle 4 is known t o have been anomalou sly long anyway (the so-called "phase 



catastrophe" in the solar cycle, Vitinskv et all Il986l ) and its decaying phase is not well 



covered by observations ( Vaauerol. 2007h. this gave rise to the suggestion of a "lost solar 



cycle" between cycles 4 and 5 ( Usoskin et al . 2001 ; Usoskin and Mursulal . l2003h . This cycle, 
however, would have been even more anomalous than cy cle 4 and despi t e inte nsive searches 



in historic data the e vidence is still not quite conclusive (jKrivova et all. 120021 : see, however, 



Usoskin et all . l2009al ). 



The issue whether the even-odd rule can go through phase jumps or not is important with 
respect to its possible origin. One plausible possibility is that the alternation is due to the super- 
position of a ste ady primordial magnetic field component on the oscillatory magnetic field generated 
by the dynamo ( Levy and Boyerl . Il982l ). In this case, any phase jump in the Gnevyshev-Ohl rule 



should imply a phase jump in Hale's polarity rules, too. Alternatively, pe rsistent even-o dd al- 
ternation may also arise in nonlinear dynamos as a period- 2 limit cycle (Durnevl . l2000h ; with 



a sto chastic forcing occasional phase jumps are also possible Icharbonneau 



2001; Charbonneau 



et al. . 120071 ). 

While we have no information on this from the 18th century phase jump, we can be certain 
that there was no such phase jump in polarities in the last two decades, even though the even-odd 
rule seems to have been broken again. It will be interesting to see when (and if) the even-odd rule 
settles in again, whether it will have done so with a phase jump or not. For instance, if cycle 25 will 
again exceed cycle 24 it would seem that no phase jump occurred and both theoretical options are 
still open. But if cycle 25 will represent a further weakening from cycle 24, followed by a stronger 
cycle 26, a phase jump will have occurred, which may exclude the primordial field origin of the 
rule if Hale's polarity rules remain unchanged. 

Let us now discuss the stronger interpretation of t he even-odd rule. In the fi rst quantitative 
study of the relative amplitudes of consecutive cycles, iGnevvshev and "Ohl (jl948h found a rather 
tight correlation between the time integrated amplitudes of even and subsequent odd cycles, while 
the correlation between odd cycles and subsequent even cycles was found to be much less strong. 
This gave rise to the notion that solar cycles come in "two-packs" as even-o dd pairs. Nagovitsyn 
et al ( 2009f ) confirmed this puzzling finding on the basis of data covering the whole period of 



telescopic observations (and renumbering cycles before 1790 in accordance with the lost cycle 
hypothesis); they also argue that cycle pair 22-23 does not deviate strongly from the even-odd 
correlation curve so it should not be considered a "real" exception to the even-odd rule. 

The fact that shortly af t er its formulation by IGnevvshev and Oh l (1948), the (strong) even-odd 
rule was used by Kopeckv ( 195dl ) to successfully predict the unusually strong cycle 19 made this 
method particularly popular for forecast purposes. However, forecasts based on the even -odd rule 
comp letely failed for cycle 23, overpredicting the amplitude by > 50% (see review bv iLi et~al. . 
20011 ) . Taken together with the implausibility of the suggested two-pack system, this shows that 
it is probably wiser to take the position that "extraordinary claims need extraordinary evidence" 
- which is yet to be provided in the case of the "strong" even-odd rule. 

Finally, i n the context of th e even-odd rule, it is also worth mentioning the three-cycle regularity 
proposed bv I Ahluwalial (|1998l ). Even though the evidence presented for the alleged triadic pattern 
is not overwhelming, this method resulted in one of the few successful predictions for the amplitude 
of cycle 23. 
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3.2.3 The Gleissberg cycle 



Besides the changes in the length of the 11-year cycle related to the amplitude-cycle length cor- 
relation, ev e n mo re significant are the variations in the period of the so-called Gleissberg cycle 
( Gleissberg! . 19391 ). This "cycle", corresponding to the 60-120 year "plate au" in Figure [5] was 
actually first noticed by Wolf, who placed it in the range 55-80 years (see iRichardl . 120041 . for a 
discussion of the history of the studies of the Gleissberg cycle). Researchers in the middle of the 
20th century characterized it as an 80-100 year variation. Figure [S] explains why so widely differing 
periods were found in different studies: the period has in fact shown a secular increase in the past 
300 years, from about 50 years in the early 18th century, to a current value exceeding 140 years. 
This incr eased length of the G leissber g evele also agrees with the results of Forgacs-Daika and 
Borkovitsl2007|). 

The detection of ~ 100 year periods in a data set of 300 years is of course always questionable, 
especially if the period is even claimed to be varying. However, the very clear, and, most impor- 
tantly, nearly linear secular trend seen in Figure [9] argues convincingly for the reality of the period 
in question. This clear appearance of the period is due to the carefully optimized choice of the 
kernel function in the time-frequency analysis , a method resulting in a so-called pseudo-Wigner 
distribution (PWD). In addition, in their study iKollath and Olahl ( 20091 ) present an extremely con- 
scientious test of the reliability of their methods, effectively proving that the most salient features 
in their PWD ar e not artefacts. (The method was subsequently also applied to stellar activity, 
Olah et all 120091 .) This is the most compelling evidence for the reality of the Gleissberg cycle yet 
presented. 



3.2.4 Supersecular cycles 

For the 210-year Suess cycle. iMcCracken and Beer (2008) present further evidence for the tempo- 
rally intermittent nature of this marked peak in the spectrum of solar proxies. The Suess cycle 
seems to have a role in regulating the recurrence rate of grand minima. Grand minima, in turn, 
only seem to occur during < 1 kiloyear intervals ("Sporer events") around the minimum of the 
~ 2400-year Hallstatt cycle. 

For further discu ssion of long t e rm va riations in solar activity we refer the reader to the reviews 
bv lBeer et all (|2006f ) and lUsoskinl (j2008f ). 



3.3 Nonlinear methods 



"...every complicated question has a simple answer which is wrong. Analysing 
a time series with a nonlinear approach is dehnitely a complicated problem. 
Simple answers have been repeatedly offered in the literature, quoting numer- 
ical values for attractor dimensions for any conceivable sy stem. 



(jHegger et all . Il999l ) 



The nonlinearities in the dynamo equations readily give rise to chaotic behaviour of the solu- 
tions. The long term behaviour of solar activity, with phenomena like grand minima and grand 
maxima, is also suggestive of a chaotic system. While chaotic systems are inherently unpredictable 
on long enough time scales, their deterministic nature does admit forecast within a limited range. 
It is therefore natural to explore this possibility from the point of view of solar cycle prediction. 



3.3.1 Attractor analysis and phase space reconstruction: the pros ... 

Assuming that the previous (M — 1) values of the sunspot number do in some way determine 
the current expected value, our problem becomes restricted to an M-dimensional phase space, the 
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dimensions being the current value and the (M — 1) previous values. With a time series of length 
N, we have N — M + 1 points fixed in the phase space, consecutive points being connected by a 
line. This phase space trajectory is a sampling of the attractor of the physical system underlying 
the solar cycle (with some random noise added to it). The attractor represents a mapping in 
phase space which maps each point into the one the system occupies in the next time step: if this 
mapping is known to a good approximation, it can be used to extend the trajectory towards the 
future. 

For the mapping to be known, M needs to be high enough to avoid self-crossings in the phase 
space trajectory (otherwise the mapping is not unique) but low enough so that the trajectory still 
yields a good sampling of the attractor. The lowest integer dimension satisfying these conditions 
is the embedding dimension D of the attractor (which may have a fractal dimension itself). 

Once the attractor has been identified, its mathematical description may be done in two ways. 
(1) Parametric fitting of th e attractor mapping in phase space. The simplest method is the 
piecewise linear fit s uggested bvlFarmer and Sidorowichl (|1987T ) and applied in several solar predic- 
tion att empts, e.g., Kurths and Ruzmaikinl (|1990i ). Using a method belonging to this group, for 
cycle 24 iKilcik et all (|2009h predict a peak amplitude of 87 to be reache d in late 2012. Alterna- 
tively, a global nonlinear fit can also be used: this is the method applied bv lSerre and Nesme-Ribesl 
( 20001 ) as the first step in their global flow reconstruction (GFR) approach. 



(2) Nonparametric fitting. The simplest nonparametric fit is to find the closest known attractor 
point to ours (in the (M — 1 )-dimension al subspace excluding the last value) and then using this 
for a prediction, as donebv ljensen1(|l993l ). (This resulted in so large random forecast errors that it 
is practically unsuitable for prediction.) Neural networks, discussed in more detail in Section [3331 
below, are a much more sophisticated nonparametric fitting device. 

(3) Indirectly, one may try to find a set of differential equations describing a system that gives 
rise to an attractor with properties similar to the observed. In this case there is no guarantee that 
the derived equations will be unique, as an alternative, completely different set may also give rise to 
a very similar attractor. This arbitrariness of the choice is not necessarily a problem from the point 
of view of prediction as it is only the mapping (the attractor structure) th at matters. Such phase 
space rec onstruction by a set o f governing equations was performed, e.g., bv lSerre and Nesme-Ribes 



(|2000h or lAguirre etaU (|2008h : "or cycle 24 the latter authors predict a peak amplitude of 65 ± 16. 



On the other hand, instead of putting up with any arbitrary set of equations correctly reproducing 
the phase space, one might make an effort to find a set with a structure reasonably similar to the 
dynamo equations so they can be given a meaningful physical interpretation. Methods following 
this latter approach will be discussed in Sections 14.41 and 14.51 



3.3.2 



the cons 



Finding the embedding dimension and the attractor structure is not a trivial task, as shown by 
the widely diverging results different researchers arrive d at. One way to fin d the correct embed- 
ding dimension is the false nearest neighbours method ( Kennel et al . 19921 ). essentially designed 
to identify self-crossing in the phase space trajectory, in which case the dimension M needs to 
be increased. But self-crossings are to some extent inevitable, due to the stochastic component 
superimposed on the deterministic skeleton of the system. 

As a result, the determination of the minimal necessary embedding dimension is usually done 
indirectly. One ind i rect m ethod fairly popular in the solar community is the approach proposed by 
ISugihara and May ( 1990f ) where the correct dimension is basically figured out on the basis of how 
successfully the model, fit to the first part of the data set, can "predict" the second part (using a 
piecewise linear mapping). 



Another widely used approach, due to lGrassberger and Procaccial (|1983l ). starts by determining 
the correlation dimension of the attractor, by simply counting how the number of neighbours in an 
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embedding space of dimension M 3> 1 increases with the distance from a point. If the attractor 
is a lower dimensional manifold in the embedding space and it is sufficiently densely sampled by 
our data then the logarithmic steepness d of this function should be constant over a considerable 
stretch of the curve: this is the correlation dimension d. Now, we can increase M gradually and 
see at what value d saturates: that value determines the attractor dimension, while the value of 
M where saturation is reached yields the embedding dimension. 

The first nonlinear time series studies of solar activity indicators suggested a time series spacing 
of 2-5 ye ars, an attractor dimensi on ~ 2-3 and an embedding dimensio n of 3-4 ( Kurths and 



Ruzmaikin, 1990l : Gizzatullina et all 1990| ). Other researchers, however, were unable to confirm 
these results, eit her reporting very different values or not finding any evidence for a low dimensional 



attractor at all ( Calvo et all 119951 : iPrice et all Il992i : ICarbonell et ~al\. Il994l : Kilci k et al 
Hanslmeier and Brajsal201C ). In particular, I would like to call attention to the paper by 



20091: 



Jensenl 



(|l993l ). which, according to ADS and WoS, has received a grand total of zero citations (!) up 
to 2010, yet it displays an exemplary no-nonsense approach to the problem of sunspot number 
prediction by nonlinear time series methods. Unlike so many other researchers, the author of that 
paper does not fool himself into believing to see a linear segment on the logarithmic correlation 
integral curve (his Figure 4); instead, he demonstrates on a simple example that the actual curve 
can be perfectly well reproduced by a simple stochastic process. 

These contradictory results obviously do not imply that the mechansim generating solar activity 
is not chaotic. For a reliable determination a long time series is desirable to ensure a sufficiently 
large number of neighbours in a phase space volume small enough compared to the global scale of 
the attractor. Solar data sets (even the cosmogenic radionuclide proxies extending over millennia 
but providing only a decadal sampling) are typically too short and sparse for this. In addition, 
clearly distinguishing between the phase space fingerprints of chaotic and stochastic processes is 
an unsolved problem of nonlinear dynamics which is not unique to solar physics. A number of 
methods have been suggested to identify chaos unambiguously in a time series but none of them 
has bee n generally ac c epted and this topic is currently a subject of ongoing research - see, e.g., the 
work of Freitas et all (|2009l ) which demonstrates that the method of "noise titration" , somewhat 
akin to the Sugihara-May algorithm, is uncapable of distinguishing superimposed coloured noise 
from intrinsically chaotic systems. 



3.3.3 



and the upshot 



Starting from the 1980s many researchers jumped on the chaos bandwagon, applying nonlinear 
time series methods designed for the study of chaotic systems to a wide variety of empirical data, 
including solar acti vity parameters. From the 1990s, however, especially after the publication of the 
influential book bv lKantz and Schreiberl (|l997h . it was increasingly realized that the applicability of 
these nonlinear algorithms does not in itself prove the predominantly chaotic nature of the system 
considered. In particular, stochastic noise superposed on a simple, regular, deterministic skeleton 
can also give rise to phase space characteristics that are hard to tell from low dimensional chaos, 
especially if strong smoothing is applied to the data. As a result, the pendulum has swung in the 
opposite direction and currently the prevailing vie w is that there is no clear cut evidence for chaos 
in solar activity data ([Panchev and Tsekovl . 120071 ) . 

One might take the position that any forecast based on attractor analysis is only as good as 
the underlying assumption of a chaotic system is: if that assumption is unverifiable from the data, 
prediction attempts are pointless. This, however, is probably a too hasty judgment. As we will 
see, the potentially most useful product of phase space reconstruction attempts is the inferences 
they allow regarding the nature of the underlying physical system (chaotic or not), even offering 
a chance to constrain the form of the dynamo equations relevant for the Sun. As discussed in 
the next section, such truncated models may be used for forecast directly, or alternatively, the 
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insight they yield into the mechanisms of the dynamo may be used to construct more sophisticated 
dynamo models. 



3.3.4 Neural networks 

Neural networks are algorithms built up from a large number of small interconnected units ( "neu- 
rons" or "threshold logic units" ) , each of which is only capable of performing a simple nonlinear 
operation on an input signal, essentially described by a step function or its generalized (rounded) 
version, a sigmoid function. To identify the optimal values of thresholds and weights parameteriz- 
ing the sigmoid functions of each neuron, an algorithm called "back propagation rule" is employed 
which minimizes (with or without human guidance) the error between the predicted and observed 
values in a process called "training" of the network. Once the network has been correctly trained, 
it is capable of further predictions. 

The point is that any arbitrary multidimensional nonlinear mapping may be approximated by 
a combination of stepfunctions to a good degree - so, as mentioned in Section 13.3.11 above, the 
neural network can be used to find the nonlinear mapping corresponding to the attractor of the 
given time series. 

More detailed introductions to the method are given bv Blais and Mertzl (l200lf) and bv Calvo 



et al. (|l995l ): the latter authors were also the first to apply a neural network for sunspot number 



prediction. Unfortunately, despite their claim of being able to "predict" (i.e. postdict) some earlier 
cycles correctly, their prediction for cycle 23 was off by a wide margin (predicted peak amplitude 
of 166 as opposed to 121 observed). One of the neural network forecasts for cvcle 24 f Maris and 
Oncical2006) was equally far o f f, pre dicting a maximum of 145 as early as December 2009, while 



another one ( Uwamahoro et al. . 20091 ) yields a more conservative value of 117.5 ± 8.5 for 2012. 
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4 Model-Based Predictions 



Progress in dynamo theory is extremely difficult, as it can be made only by un- 
derstanding the interaction of turbulent plasma motions with magnetic fields. 
Indeed, the extreme conditions within the solar interior make this a formidable 
task. ... Any predictions made with such models should be treated with ex- 
treme caution (or perhaps disregarded), as they lack solid physical underpin- 
nings. 



([Tobias et ad I2006T ) 



While attempts to predict future solar cycles on the basis of the empirical sunspot number record 
have a century-old history, predictions based on physical models of solar activity only started a 
few years ago. The background of this new trend is, however, not some significant improvement 
in our understanding of the solar dynamo. Rather, it is the availability of increasingly fast new 
computers that made it possible to fine-tune the parameters of certain dynamo models to reproduce 
the available sunspot record to a good degree of accuracy and to apply data assimilation methods 
(such as those used in terrestrial weather prediction) to these models. This is not without perils. 
On the one hand, the capability of multiparametric models to fit a multitude of observational data 
does not prove the conceptual correctness of the underlying model. On the other hand, in chaotic 
or stochastic systems such as the solar dynamo, fitting a model to existing data will not lead to a 
good prediction beyond a certain time span, the extent of which can only be objectively assessed by 
"postdiction" tests, i.e. checking the models predictive skill by trying to "predict" previous solar 
cycles and comparing those predictions to available data. Apparently successful postdiction tests 
have led s ome groups to claim a breakthrough in solar cycle prediction owing to the model-based 
approach (jDikpati and Gilmanl . 2006; iKitiashvili and Kosovichevl . l2008l ). Yet, as we will see in the 



following discussion, a closer inspection of these claims raises many questions regarding the role 
that the reliance on a particular physical dynamo model plays in the success of their predictions. 



4.1 The solar dynamo: a brief summary of current models 

Ex tensive summari es of the current sta n ding of solar dynamo theor y are given in the re views 
bv lPetrovavl (|2000l ). lOssendriiverl (|2003l ). ICharbonneaul (j2005l ). and fsolanki et all (|2006l ). As 
explained in detail in those reviews, all current models are based on the mean-field theory approach 
wherein a coupled system of partial differential equations governs the evolution of the toroidal and 
poloidal components of the large-scale magnetic field. The large-scale field is assumed to be axially 
symmetric in practically all current models. In some nonlinear models the averaged equation of 
motion, governing large-scale flows is also coupled into the system. 

In the simplest case of homogeneous and isotropic turbulence, where the scale I of turbulence 
is small compared to the scale L of the mean variables (scale separation hypothesis), the dynamo 
equations have the form 



dB 

~dt 



= V x (U x B + aB) - V x (?y T x VB) 



(13) 



Here B and U are the large-scale mean magnetic field and flow speed, respectively; t]t is the 
magnetic diffusivity (dominated by the turbulent contribution for the highly conductive solar 
plasma), while a is a parameter related to the non-mirror symmetric character of the magnetized 
plasma flow. 

In the case of axial symmetry the mean flow U may be split into a meridional circulation U c 
and a differential rotation characterized by the angular velocity profile 9): 



U = U r 



r sin 9 fio e<i 
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where r, 9, </> are spherical coordinates and is the azimuthal unit vector. Now introducing the 
shear 

CI = r sin 9 Vf2 , £1 = — sgn — — ■ H , 

dr 

assuming a -C f2£ and ignoring spatial derivatives of a and r/T, Equation (|13p simplifies to the 
pair 

— =a J B-(U c -V)A-(V-U c )A + ry T V 2 A (14) 

— = — - (U c • V)B - (V • U C )B + VT V 2 B, (15) 

where B and A are the toroidal (azimuthal) components of the magnetic field and of the vector 
potential, respectively, and is to be evaluated in the direction 90° clockwards of il (along the 
isorotation surface) in the meridional plane. These are the classic aQ dynamo equations, including 
a meridional flow. 

In the more mainstream solar dynamo models the strong toroidal field is now generally thought 
to reside near the bottom of the solar c onvective zone. Indeed, it is known that a variety of flux 
transport mechanisms such as pumping ( Petrovav , 1994f ) remove magnetic flux from the solar con- 
vective zone on a timescale short compared to the solar cycle. Following earlier simpler numerical 
experiments, recent MHD numerical simulations have indeed demonstrated this pumping of large 
scale magnetic flux fro m the convective zone i nto the tachocline below, where it forms strong 



coherent toroidal fields (jBrowning et all I2006T ) . As this layer is also where rotational shear is 



maximal, it is plausible that the strong toroidal fields are not just stored but also generated here, 
by the winding up of poloidal field. The two main groups of dynamo models, interface dynamos 
and flux transport dynamos, differ mainly in their assumptions about the site and mechanism of 
the a-effect responsible for the generation of a new poloidal field from the toroidal field. 

In interface dynamos a is assumed to be concentrated near the bottom of the convective zone, 
in a region adjacent to the tachocline, so that the dynamo operates as a wave propagating along the 
interface between these two layers. While these models may be roughly consistent and convincing 
from the physical point of view, they have only had limited success in reproducing the observed 
characteristics of the solar cycle, such as the butterfly diagram. 

Flux transport dynamos, in contrast, rely on the Babcock-Leighton mechanism for a, arising 
due to the action of the Coriolis force on emerging flux loops, and they assume that the corre- 
sponding a-effect is concentrated near the surface. They keep this surface layer incommunicado 
with the tachocline by introducing some arbitrary unphysical assumptions (such as very low diffu- 
sivities in the bulk of the convective zone) . The poloidal fields generated by this surface a-effect 
are then advected to the poles and there down to the tachocline by the meridional circulation - 
which, accordingly, has key importance in these models. The equatorward deep return flow of 
the meridional circulation is assumed to have a significant overlap with the tachocline (another 
controversial point), and it keeps transporting the toroidal field generated by the rotational shear 
towards the equator. By the time it reaches lower latitudes, it is amplified sufficiently for the 
flux emergence process to start, resulting in the formation of active regions and, as a result of 
the Babcock-Leighton mechanism, in the reconstruction of a poloidal field near the surface with 
a polarity opposed to that in the previous 11-year cycle. While flux transport models may be 
questionable from the point of view of their physical consistency, they can be readily fine-tuned to 
reproduce the observed butterfly diagram quite well. 

It should be noted that while the terms "interface dynamo" and "flux transport dynamo" are 
now very widely used to describe the two main approaches, the more generic terms "advection- 
dominated" and "diffusion-dominated" would be preferable in several respects. This classification 
allows for a continuous spectrum of models depending on the numerical ratio of advective and 
diffusive timescales (for communication between surface and tachocline). In addition, even at 
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the two extremes, classic interface dynamos and circulation-driven dynamos are just particular 
examples of advection or diffusion dominated systems with different geometrical structures. 



4.2 Is model-based cycle prediction feasible? 

As it can be seen even from the very brief and sketchy presentation given above, all current solar 
dynamo models are based on a number of quite arbitrary assumptions and depend on a number of 
free parame ters, the functional form a nd amplitude of which is far from being well constrained. For 
this reason. iBushby and Tobias! (|200 7f ) rightfully say that all current solar dynamo models are only 



of "an illustrative nature" . This would suggest that as far as solar cycle prediction is concerned, 
the best we should expect from dynamo models is also an "illustrative" reproduction of a series 
of solar cycles with the same kind of long-term varia tions (qualitatively and, in the statistical 



sense, quantitatively) as seen in solar data. Indeed, iBushbv and Tobias! (|2007l ) demonstrated 
that even a minuscule stochastic variation in the parameters of a particular flux transport model 
can lead to large, unpredictable variations in the cycle amplitudes. And even in the absence of 
stochastic effects, the chaotic nature of nonlinear dynamo solutions seriously limits the possibilities 
of prediction, as the authors find in a particular interface dynamo model: even if the very same 
model is used to reproduce the results of one particular run, the impossibility of setting initial 
conditions exactly representing the system implies that predictions are impossible even for the 
next cycle. Somewhat better results are achieved by an alternative method, based on the phase 
space reconstruction of the attractor of the nonlinear system - this is, however, a purely empirical 
time series analysis technique for which no knowledge of the detailed underlying physics is needed. 
(Cf. Section 1531 above.) 

Despite these very legitimate doubts regarding the feasibility of model-based prediction of solar 
cycles, in recent years several groups have claimed to be able to predict solar cycle 24 on the basis 
of dynamo models with a high confidence. So let us consider these claims. 

4.3 Explicit models 

The current buzz in the field of model-based solar cycle prediction was started by the work of the 
solar dynamo group in Boulder ( Dikpati et al. . 2006k Dikpati and Gilman . 20061 ). Their model is a 



flux transport dynamo, advection-dominated to the extreme. The strong suppression of diffusive 
effects is assured by the very low value (less than 20km 2 /s) assumed for the turbulent magnetic 
diffusivity in the bulk of the convective zone. As a result, the poloidal fields generated near 
the surface by the Babcock-Leighton mechanism are only transported to the tachocline on the 
very long, decadal time scale of meridional circulation. The strong toroidal flux residing in the 
low-latitude tachocline, producing solar activity in a given cycle is thus the product of the shear 
amplification of poloidal fields formed near the surface about 2-3 solar cycles earlier, i.e. the 
model has a "memory" extending to several cycles. The mechanism responsible for cycle-to-cycle 
variation is assumed to be the stochastic nature of the flux emergence process. In order to represent 
this variability realistically, the model drops the surface a-term completely (a separate, smaller 
a term is retained in the tachocline); instead, the generation of poloidal field near the surface is 
represented by a source term, the amplitude of which is based on the sunspot record, while its 
detailed functional form remains fixed. 

Dikpati and Gilman (|2006l ) find that, starting off their calculation by fixing the source term 



amplitudes of sunspot cycles 12 to 15, they can predict the amplitudes of each subsequent cycle 
with a reasonable accuracy, provided that the relation between the relative sunspot numbers and 
the toroidal flux in the tachocline is linear, and that the observed amplitudes of all previous cycles 
are incorporated in the source term for the prediction of any given cycle. For the upcoming cycle 24 
the model predicts peak smoothed annual relative sunspot numbers of 150 or more. Elaborating 
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on their model, they proceeded to apply it separately to the northern and southern hemispheres, 
to find t hat the model can a lso be used to correctly forecast the hemispheric asymmetry of solar 
activity (jDikpati et adl2007t ). 



The extraordinary claims of this pioneering research have prompted a hot debate in the dynamo 
community. Besides the more general, fundamental doubt regarding the feasibility of model-based 
predictions (see Section 3.2 above), more technical concerns arose, to be discussed below. 

Another flux transport dynamo code, the Surya code, originally developed by A. Choudhuri 
and coworkers in Bangalore, has also been utilized for prediction purposes. The crucial difference 
between the two models is in the value of the turbulent diffusivity assumed in the convective 
zone: in the Bangalore model this value is 240 km 2 /s, 1-2 orders of magnitude hig her than in the 



Boulder model, and within the physically plausible range (jChatteriee et all 120041 ) . As a result of 
the shorter diffusive timescale, the model has a shorter memory, not exceeding one solar cycle. As 
a consequence of this relatively rapid diffusive communication between surface and tachocline, the 
poloidal fields forming near the surface at low latitudes due to the Babcock-Leighton mechanism 
diffuse down to the tachocline in about the same time as they reach the poles due to the advection 
by the meridional circulation. In these models, then, polar magnetic fields are not a true physical 
precursor of the low-latitude toroidal flux, and their cor relation is just due to their common source . 



In the version of the code adapted for cycle prediction (jChoudhuri et al. , 2007 ; Jiang et al. L l2007h . 



the "surface" poloidal field (i.e. the poloidal field throughout the outer half of the convection 
zone) is rescaled at each minimum by a factor reflecting the observed amplitude of the Sun's dipole 
field. The model shows reasonable predictive skill fo r the last three cycle s for w hich data are 
available, and can even tackle hemispheric asymmetry ( Goel and Choudhuri . 20091 ). For cycle 24, 
the predicted amplitude is 30-35% lower than cycle 23. 

4.4 Truncated models 

The "illustrative" nature of solar dynamo models is nowhere more clearly on display than in 
truncated or reduced models where some or all of the detailed spatial structure of the system is 
completely disregarded, and only temporal variations are explicitly considered. This is sometimes 
rationalised as a truncation or spatial integration of the equations of a more realistic inhomogenous 
system; in other cases, no such rationalisation is provided, representing the solar dynamo by an 
infinite, homogeneous or periodic turbulent medium where the amplitude of the periodic large-scale 
magnetic field varies with time only. 

In the present subsection we deal with models that do keep one spatial variable (typically, the 
latitude), so growing wave solutions are still possible - these models, then, are still dynamos even 
though their spatial structure is not in a good correspondence with that of th e solar dynamo . 

This approach in fact goes back to the classic migratory dynamo model of iParkerl (|1955f ) who 
radially truncated (i.e., integrated) his equations to simplify the problem. Parker seems to have 
been the first to employ a heuristic relaxation term of the form —B r /Td in the poloidal field 
equation to represent the effect of radial diffusion; here, = d 2 /r)T is the diffusi ve timescale 
across the thickness d of the conv ective zone. His model was recently generalized by M oss et al 



(2008) and lUsoskin et ah (|2009bl ) ;o the case when the a-effect includes an additive stochastic 



noise, and nonlinear saturation of the dynamo is achieved by a-quenching. These authors do not 
make an attempt to predict solar activity with their model but they can reasonably well reproduce 
some features of the very long term solar activity record, as seen from cosmogenic isotope studies. 

Another radially truncated model, this time formulated in a Cartesian svstcm. is that of Kiti- 
ashvili and Kosovichev (|2009l ). In this model stochastic effects are not considered, and, in addition 
to using an a-quenc hing recipe, further np nlinearity is introduced by coupling in the Kleeorin- 
Ruzmaikin equation (jZeldovich et al. , 1984h governing the evolution of magnetic helicity, which in 



the hydromagnetic case contributes to a. Converting the toroidal field strength to relative sunspot 
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number using the Bracewell transform, eq. ([3]), the solutions reproduce the asymmetric profile of 
the sunspot number cycle. For sufficiently high dynamo numbers the solutions become chaotic, 
cycle amplitudes show an irregular variation. Cycle amplitudes and minimum-maximum time 
delays are found to be related in a way reminiscent of the Waldm eier relation. 

Building on these results, Kitiashvili and Kosovichev ( 20081 ) attempt to predict solar cycles 
using a data assimilation method. The approach used is the so-called Ensemble Kalman Filter 
method. Applying the model for a "postdiction" of the last 8 solar cycles yielded astonishingly 
good results, considering the truncated and arbitrary nature of the model and the fundamental 
obstacles in the way of reliable prediction discussed above. While the presently available brief 
preliminary publication leaves several details of the method unclear, the question may arise whether 
the actual physics of the model considered has any significant role in this prediction, or we are 
dealing with something like the phase space reconstruction approach discussed in Section [3.31 above 
where basically any model with an attractor that looks reasonably similar to that of the actual 
solar dynamo would do. Either way, the method is remarkable, and the prediction for cycle 24 of a 
maximal smoothed annual sunspot number of 80, to be reached in 2013, will be worth comparing 
to the actual value. 

In or der to understand the origin of the predictive sk ill of the Boulde r model. Cameron and 
Schiisslerl2007|) studied a radially truncated version of the model, wherein only the equation for 
the radial field component is solved as a function of time and latitude. The equation includes a 
source term similar to the one used in the Boulder model. As the toroidal flux does not figure 
in this simple model, the authors use the transequatorial flux $ as a proxy, arguing that this 
may be more closely linked to the amplitude of the toroidal field in the upcoming cycle than 
the polar field. They find that $ indeed correlates quite well (correlation coefficients r ~ 0.8-0.9, 
depending on model details) with the amplitude of the next cycle, as long as the form of the latitude 
dependence of the source term is prescribed and only its amplitude is modulated with the observed 
sunspot number series ("idealized model"). But surprisingly, the predictive skill of the model is 
completely lost if the prescribed form of the source function is dro pped and the actually observed 
latitude distribution of sunspots is used instead ("realistic model"). ICameron and Schiissler ( 2007 ) 
interpret this by pointing out that $ is mainly determined by the amount of very low latitude flux 
emergence, which in turn occurs mainly in the last few years of the cycle in the idealized model, 
while it has a wider temporal distribution in the realistic model. The conclusion is that the root 
of the apparently good predictive skill of the truncated model (and, by inference, of the Boulder 
model it is purported to represent) is actually just the good empirical correlation between late- 
phase activity and the amplitude of the next cycle, discussed in Section [2 . 1 1 above . This correlation 
is implicitly "imported" into the idealized flux transport model by assuming that the late-phase 
activity is concentrated at low latitudes, and therefore gives rise to cross-equatorial flux which 
then serves as a seed for the toroidal field in the next cycle. So if Cameron and Schiisslerl (|2007l ) 
are correct, the predictive skill of the Boulder model is due to an empirical precursor and is thus 
ultimately explained by the good old Waldmeier effect (cf. Sect i on|1.3.3j) 

The fact that the truncated model of lCameron and Schiisslerl (|2007l ) is not identical to the Boul- 
der model obviously leaves room for doubt regarding this conclusion. In particular, the effective 
diffusivity represented by the sink term in the truncated model is ~ 100 km 2 /s, significantly higher 
than in the Boulder model; consequently, the truncated model will have a more limited memory, cf. 
Yeates et all ( 20081 ). The argument that the cross-equatorial flux is a valid proxy of the amplitude 



of the next cycle may be correct in such a short-memory model with no radial structure, but it is 
dubious whether it remains valid for flux transport models i n general. In a n attem pt to appreciate 
the importance of the cross-equatorial flux in their model, iDikpati et all (|2008al ) find that while 
this flux does indeed correlate fairly well (r = 0.76) with the next cycle amplitude, the toroidal 
flux is a much better predictor (r = 0.96). At first sight this seems to make it unlikely that the 
former can explain the latter; however, part of the difference in the predictive skill may be due to 
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the fact that $ shows much more short-term variability than the t oroidal flux. 

In any case , the obvious wa y to address the concerns raised by ICameron and Schussler ( 2007 ) 



and further bv lSchiisslerl ([20071 ) in relation to the Boulder model would be to run that model with 
a modified source function incorporating the realistic latitudinal distribution of sunspots in each 
cycle. The results of such a test are not yet available at the time of writing this review. 

4.5 The Sun as an oscillator 

An even more radical simplification of the solar dynamo problem ignores any spatial dependence in 
the solutions completely, concentrating on the time dependence only. Spatial derivatives appearing 
in Equations (fl"4|) and (|15[) are estimated as V ~ 1/L and the resulting terms U c /L and t]t/L 2 as 
1/r where r is a characteristic time scale. This results in the pair 

A = aB-A/r (16) 

B={n/L)A-B/r (17) 

which can be combined to yield 

.. D - 1 2 ■ 

B= — —B--B (18) 
r z r 

where D = aflr 2 /L is the dynamo number. For D < 1, Equation (p~8]> clearly describes a damped 
linear oscillator. For D > 1, solutions have a non-oscillatory character. The system described by 
Equation (fT8|) . then, is not only not a true dynamo (missing the spatial dependence), but it does 
not even display growing oscillatory solutions that would be the closest counterpart of dynamo- 
like behaviour in such a system. Nevertheless, there are a number of ways to extend the oscillator 
model to allow for persistent oscillatory solutions, i.e., to turn it into a relaxation oscillator. 

(1) The most straightforward approach is to add a forcing term + sin(o;ot) to the r.h.s. of 
Equations (fl~8|) . Damping would cause the system to relax to t h e driy ing period 2tt/u>o if there 
were no stochastic disturbances to this equilibrium. iHiremathl (120061) fitted the parameters of 



the fo rced and dampe d oscillator model to each observed solar cycle individually; then in a later 
work (|Hiremathl . [2008h he applied linear regression to the resulting series to provide a forecast (see 



Section 13.11 above) . 

(2) Another trick is to account for the ir/2 phase difference between poloidal and toroidal field 
components in a dynamo wave by introducing a phase factor i into the first term on the r.h.s. of 
Equation (IT7|) . This can also be given a more formal derivation as equations of this form result 
from the substitution of solutions of the form A oc e lkx , B cx e l ( kx +^/ 2 ) into the ID dynamo 
equations. This route, combined with a nonlinearity due to magnetic modulation of differential 
rotation described by a coupled third equation, was taken by IWeiss et all (jl984h . Their model 



displa yed chaotic behaviour with intermittent episodes of low activity similar to grand minima. 
(31 IWilmot -Smith et all (120061) showed that another case where dynamo-like behaviour can be 



found in an equation like (|18l) is if the missing effects of finite communication time between parts 
of a spatially extended system are reintroduced by using a time delay At, evaluating the first term 
on the r.h.s. at time t — At to get the value for the l.h.s. at time t. 

(4) Yet another possibility is to introduce a nonlinearity into the model by assuming D = 
Do[l — f(B)] where f(B = 0) = and / > everywhere. (Note that any arbitrary form of a- 
or f2-quenching can be cast in the above form by series expansion.) The governing equation then 
becomes one of a nonlinear oscillator: 

B = ^±B --B- ^——Bf{B) (19) 

In the most commonly assumed quenching mechanisms the leading term in f(B) is quadratic 



In the most commonly assumed quenching mechanisms the leading term m f{±>) is quadratic; 
in this case Equation (TT9l describes a Duffing oscillator (jKanamarul . 2008f ). For large positive 
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dynamo numbers, Dq 3> 1, then, the large nonlinear term dominates for high values of B, its 
negative sign imposing oscillatory behaviour; yet the origin is a repeller so the oscillation will 
never be damped out. The Duffing oscillator was first considered in the solar context by Palus and 



Novotna ( 19991). Under ce r tain condition s on the parameters, it c an be reduced to a van der Pol 
oscillator ( Adomian . 19891 : Mininni et at . 2002 : Kanamaru . 2007 ): 



(20) 



with fi > 0. From this form it is evident that the problem is equivalent to that of an oscillator with 
a damping that increases with amplitude; in fact, for small amplitudes the damping is negative, 
i.e., the oscillation is self-excited. 

These simple nonlinear oscillators were among the first physical systems where chaotic be- 
haviour was detected (when a periodic forcing was added). Yet, curiously, they first emerged in 
the solar context precisely as an alternative to chaotic behaviour. Consider ing the mapping of 
the solar cycle in the differential phase space {B , dB / dt} ^MmiImi et all ( 2000l ) got the impression 
that, rather than showing signs of a strange attractor. the SSN series is adequately modelled 
by a van der Pol o scillat or with stochastic fluctuations. This concept was further developed by 
lLopes and Passosl ( 20091 ) who fitted the parameters of the oscillator to each individual sunspot 
cycle. The parameter fj, is related to the meridional flow speed and the fit indicates that a slower 
meridional flow may have been responsible for the Dalton minimum. This was also corroborated 
in an explicit dynamo model (the Surya code) — however, as we discussed in Sect. 12.41 this result 
of flux transport dynamo models is spurious and the actual effect of a slower meridional flow is 
likely to be opposite to that suggested by the van der Pol oscilla tor model. 

In an alternative approach to the problem, Nagovits'vnl ( 19971 ) attempted to constrain the prop- 
erties of the solar oscillator from its amplitude-frequency diagram, suggesting a Duffing oscillator 
driven at two secular periods. While his empirical reconstruction of the amplitude-frequency plot 
may be subject to many uncertainties, the basic idea is certainly noteworthy. 

In summary: despite its simplicity, the oscillator representation of the solar cycle is a relatively 
new development in dynamo theory, and its obvious potential for forecasting purposes has barely 
been exploited. 
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5 Summary Evaluation 



The perform ance of various forecast methods in cycles 21-23 was discussed bv lLi et all (|200lh and 
Kand (l200lh . 



Precursor methods stand out with their internally consistent forecasts for these cycles which for 
cycles 21 and 22 proved to be correct. For cycle 23 these methods were still internally consistent 
in their prediction, mostly scattering in a narrow range between 150 and 170; however, the cycle 
amplitude proved to be considerabl y lower (i? max = 121). It should be noted, however, that one 
precursor based prediction, that of ISchatten et~al\ (|l996f ) was significantly lower t han the res t 



(138 ± 30) and within 0.6a of the actual value. Indeed the method of Schatten et al. (|l987l : ll996T ). 
has consistently proven its skill in all cycles. As discussed in Sect. 12.21 this method is essentially 
based on the polar magnetic field strength as precursor. 

Extrapolation methods as a whole have shown a much less impressive performance. Overall, 
the statistical distribution of maximum amplitude values predicted by "real" forecasts made using 
these methods (i.e. forecasts made at or before the minimum epoch) for any given cycle does not 
seem to significantly differ from the long term climatological average of the solar cycle quoted in 
Section [L3l above (100 ±35). It would of course be a hasty judgement to dismiss each of the widely 
differing individual approaches comprised in this class simply due to the poor overall performance 
of the group. In particular, some novel methods suggested in the last 20 years, such as SSA or 
neural networks have hardly had a chance to debut, so their further performance will be worth 
monitoring in upcoming cycles. 

One group of extrapolation methods that stands apart from the rest are those based on the even- 
odd rule. These methods enjoyed a relatively high prestige until cycle 23, when they coherently 
predicted a peak amplitude around 200, i.e., ~ 70% higher than the actual peak. This can only 
be qualified as a miserable failure, independently of the debate as to whether cycle 23 is truly at 
odds with the even-odd rule or not. 

In this context it may be worth noting that the double peaked character and long duration of 
cycle 23 implies that its integrated amplitude (sum of annual sunspot numbers during the cycle) 
is much less below that of cycle 22 than the peak amplitude alone would indicate. This suggests 
that forecasts of the integrated amplitude (rarely attempted) could be more robust than forecasts 
of the peak. Nevertheless, one has to live with the fact that for most practical applications (space 
weather) it is the peak amplitude that matters most, so this is where the interest of forecasters is 
naturally focused. 

Finally, model based methods are a new development that have had no occasion yet to prove their 
skill. As discussed above, current dynamo models do not seem to be at a stage of development 
where such forecasts could be attempted w i th any confidence, especially before the time of the 



minimum. (The method of lChoudhuri et aU , 120071 using polar fields as input near the minimum, 



would seem to be akin to a version of the polar field based precursor method with some extra 
machinery built into it.) The claimed good prediction skills of models based on data assimilation 
will need to be tested in future cycles and the roots of their apparent success need to be understood. 

Table Q] presents a collect i on of forecasts for the amplitude of cycle 24, without claiming com- 
pleteness. (See, e.g., Pesnell . 120081 . for a more exhaustive list.) The objective was to include one 



or two representative forecasts from each category. 

The incipient cycle 24 may be a milestone for solar cycle forecasting. Current evidence indicates 
that we are at the end of the Modern Maximum when the Sun is about to switch to a state of less 
intense long term activity. The appearance of a number of novel prediction methods, in particular 
the model based approach, as well as the unusually large discrepancy between forecasts based on 
the precursor approach imply that, whichever course solar activity will take in the coming years, 
we have a lot to learn from the experience. 
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Table 1: A selection of forecasts for cycle 24. 



Category Peak amplitude Link Reference 



Precursor methods 








Minimax 


80 ± 25 


Eg. urn 


Brown (1976): Braisa et al. (2009)* 


Minimax3 


69 ± 15 


Eq.HTI 


Cameron and Schusslcr (2007) 


T. 1 -G 1 J 

Polar held 


75 ± 8 


Sect. 12.21 


Svalgaard et al. (2005) 


Polar ncld 


80 ± 30 


Sect. 12.21 


O 1 -L-i- /<1AAr\ 

Schattcn (2005) 


Geomagnetic (Fcynman) 


lot) 


beet. Iz.jI 


Matnawav ana Wilson (zUUbj 


Geomagnetic (Ohl) 


93 ±20 


Sect. 1231 


Bhatt et al. (2009) 


Geomagnetic (Ohl) 


101 ±5 


Sect. |3.2.2| 


Ahluwalia and Yebuhav C2009) 


Geomagnetic (interpl.) 


97 ±25 


Sect. 1231 


Wang and Sheelev ^OOg) 


Field reversal 


94 ± 14 


Eq.DJ 


Tlatov (2009)* 


Extrapolation methods 








Linear regression 


90 ±27 


Sect. O 


Braisa et al. (2009) 


Linear regression 


110 ± 10 


Sect. O 


Hiremath (2008) 


Spectral (MEM) 


90 ± 11 


Sect. OLD 


Kane (2007) 


Spectral (SSA) 


117 


Sect. E2 


Loskutov et al. (2001) 


Spectral (SSA) 


106 


Sect. I3~2l 


Kuzanvan et al. (2008) 


Attractor analysis 


87 


Sect.|3.3.1| 


Kilcik et al. (2009) 


Attractor analysis 


65 ± 16 


Sect. |3.3.1| 


Auuirre et al. (2008) 


Attractor analysis 


145 ±7 


Sect. |3.3.1| 


Crosson and Binder (2009) 


Neural network 


145 


Sect. |3.3.4| 


Maris and Oncica (2006) 


Neural network 


117.5 ±8.5 


Sect. |3.3.4| 


Uwamahoro et al. (2009) 


Model based methods 








Explicit models 


167 ± 12 


Sect. [O] 


Dikpati and Gilman (2006) 


Explicit models 


- 80 


Sect. [L3| 


Choudhuri et al. (2007) 


Explicit models 


- 85 


Sect. EM 


Jiane et al. (2007) 


Truncated models 


- 80 


Sect.SU 


Kitiashvili and Kosovichev (2008) 



References marked with * are to the basic principle used in the given prediction method while the actual numerical 
evaluation for cycle 24 was done by the author. The application for forecast purposes does not necessarily reflect 
the original intention of the basic principle, as laid out in the cited publications. 
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6 Epilogue 



Throughout the ages, mankind felt and tried to answer the urge to predict events to come. Omens 
were carefully collected and categorized on Mesopotamian clay tablets; omen-based prediction 
was developed into an industry in the form of hepatoscopy (analyzing the shape of the liver of 
a sacrificed animal) and, in later Roman times, of auspicium (watching the flight of the birds). 
Ancient Greeks often turned to oracles like the Pythia of Delphi. By the late Antiquity, the 
astrological world view was widespread throughout the civilized world, implying that cosmic and 
terrestrial events were subject to cosmic cycles governed by a variety of (planetary) periods. 

Today we tend to smile at these "superstitious" early attempts. Yet, ironically, many of the 
"advanced" methods we have for the prediction of solar activity are based on principles that 
hardly differ from those listed above: just substitute "precursor" for "omen" , "neural network" for 
"oracle" or "harmonic analysis" for "cosmic cycles"... 

But in parallel with the often naive phenomenological or empirical prediction attempts, already 
in the Hellenistic world, a handful of enlightened scientists started the development of physical mod- 
els, based on logic a nd experienc e, that would lead to the advanced predictive skills of many models 
of modern science ( Russol . 2004). Extending the analogy, we can see that the real importance of 
the recent debut of model-based solar cycles predictions is not their still dubious success rate but 
the conceptual leap they represent. 

Despite the rather poor overall performance of solar cycle prediction attempts, the extensive 
efforts invested in this endeavour were not in vain as they have contributed and keep contributing 
to a better understanding of the physical processes governing the solar cycle and to constraining 
the dynamo. 
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